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ABSTRACT 

(Distribution  Limitation  Statement  B) 


^Analysis  of  data  from  published  cratering  experiments  shows 
the  effect  of  soil  and  rock  properties  on  the  apparent  dimensions 
of  explosion  produced  craters.  ;>More  than  200  cratering  tests 
/clnd  related  material  properties  were  cataloged.  The  data  con¬ 
sisted  of  10  nuclear  events  whose  yields  varied  from  0.42  to  100 
'kilotons  and  about  200  high  explosive  events  whose  yields  varied 
from  1  to  1  million  pounds  of  TNT.  The  different  test  sites 
included  materials  for  which  the  density  ranged  from  60  to  170 
pounds/cubic  foot. 

'—^-By  regression  analysis,,  using  bell  shaped  curves,,  prediction 
formulas  were  developed  for  the  apparent  crater  radius,  depth, 
and  volume  as  a  function  of  charge  weight  and  depth  of  burst  for 
eight  different  types  of  materials.  The  bell  curves  were 
normalized  using  material  properties  and  prediction  equations 
were  generated  using  all  the  data.  These  general  equations  were 
then  studied  to  determine  the  specific  effects  of  the  material 
properties  on  resultant  apparent  crater  dimensions.  _ _ _ _ „ 


C-Meterial  properties  are  highly  important  in  determining  the 
size  of  explosion-produced  craters^and  some  of  the  more  important 
properties  are  unit  weight,  degree  of  saturation,  shearing  res'i's - ~ 
tance  and  seismic  velocity.  Previous  investigators  have  been 
somewhat  negligent  in  measuring  material  properties  for  past 
cratering  experiments  and  no  real  data  analysis  can  be  made  until 
the  variables  are  either  controlled  or  measured.  Material 
properties  which  should  be  measured  for  future  tests  should  at 
least  include  the  above  properties  and  if  possible  the  material's 
energy  dissipation  and  bulking  characteristics.  Better  yet  3 
reasonably  simple  theory  of  cratering  is  needed  which  will  better 
define  the  material  properties  governing  cratering  mechanics. 
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The  following  symbols  are  used  in  this  report: 

A.  =  coefficient  of  bell  curve  in  coded  form 

B.  =  coefficient  of  bell  curve; 
i 

C.  =  general  coefficient; 

D  =  crater  depth; 

D  =  maximum  depth  of  the  apparent  crater; 
Dob  =  depth  of  burst; 

Ev  =  vaporization  energy; 

G$  =  grain  bulk  specific  gravity; 

Ha-|  =  apparent  crater  lip  crest  height; 

La  =  linear  apparent  crater  dimension 

(radius,  depthjor  cube  root  of  volume); 
M  =  general  material  property; 

R  =  crater  radius; 

R  =  radius  of  apparent  crater; 

cl 

r  ^  =  multiple  correlation  coefficient 
m 

S  =  degree  of  saturation; 

SGZ  =  surface  ground  zero; 

TNT  =  the  high  explosive,  trinitrotluene; 

V  =  crater  volume; 

Vg  =  volume  of  apparent  crater; 

W  =  weight  of  explosive; 

xi 


NOTATION  (Continued 

ZP  =  zero  point-effective  center  of  explosion  energy; 
c  =  seismic  velocity; 
exp  =  exponential  (e); 
g  =  acceleration  due  to  gravity; 

In  =  natural  logarithm; 

m  =  the  divisor  portion  of  the  scaling  exponent,  1/m 
tan  <j>  =  tangent  of  the  angle  of  shearing  resistance; 
x,y  =  rectangular  Cartesian  coordinates; 
y  =  total  unit  weight; 

Yd  =  dry  unit  weight; 
p  =  mass  unit  weight. 


SECTION  I 


# 


INTRODUCTION 

Over  the  past  several  years,  because  of  intensified  studies 
of  possible  engineering  applications  of  nuclear  energy,  increased 
attention  has  been  devoted  to  the  problem  of  cratering  by  explo¬ 
sives.  One  of  the  most  obvious  peaceful  applications  of  nuclear 
explosives  is  that  of  earth  excavation,  as  might  be  considered 
for  the  construction  of  harbors,  dams,  or  canals.  Prediction 
and  design  for  survival  of  silo-launched  missile  systems  has 
also  added  urgency  to  these  studies. 

Although  most  investigators  have  recognized  that  the  size 
of  a  crater  obtained  from  an  explosive  charge  depends  upon  the 
media  in  which  it  is  detonated,  properties  describing  the  media 
which  relate  to  the  dimensions  of  this  crater  are  somewhat 
obscure.  The  primary  purpose  of  this  investigation,  then,  was 
to  determine  which  engineering  properties  normally  measured  for 
earth  and  rock  materials  could  be  related  to  the  size  of  a 
crater  created  by  an  explosive  charge. 

A  large  number  of  cratering  experiments  have  been  conducted 
in  various  media  (15,60,87).  These  experiments  primarily 
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provided  data  on  the  effect  of  explosive  energy  and  depth  of 
burst  on  crater  dimensions.  Although  the  experiments  number 
in  the  thousands  and  although  engineering  material  properties 
were  measured  for  a  good  percentage  of  these  experiments,  very 
little  has  been  reported  which  relate  these  material  properties 
to  the  final  crater  geometry.  Previous  investigators  seem  to 
have  been  of  two  minds:  (1)  Because  material  properties  vary 
greatly  within  one  media  in  one  location  and  because  accurate 
measurement  of  these  properties  is  often  difficult,  the  best 
approaches  are  to  ignore  completely  the  material  properties  or 
to  over  simplify  and  let  the  material  be  described  by  one 
constant  in  a  particular  prediction  relationship;  and  (2)  be¬ 
cause  the  cratering  process  is  so  complicated,  extensive 
measurement  of  material  properties  both  in  the  field  and  in 
the  lab  are  required  to  describe  the  media  to  be  subjected  to 
an  explosive  charge.  It  would  appear  that  the  better  solution 
lies  somewhere  between  these  two  extremes. 


SECTION  II 


GENERAL  CONCEPTS  AND  TERMINOLOGY 

This  sectioq  of  the  report  presents  basic  concepts  concern¬ 
ing  the  parameters  and  mechanisms  involved  in  the  explosive 
formulation  of  craters.  The  effect  of  explosive  weight,  depth 
of  burst  and  type  of  material  are  discussed.  Terminology 
generally  associated  with  the  cratering  process  is  also  presented. 

General  Description  of  the  Cratering  Process  (63).  When  an 
explosion  occurs  at  or  near  the  surface  of  a  soil  or  rock-like 
material,  a  crater  is  formed  (see  Fig.  1).  The  size  of  this 
crater  depends  on  at  least  four  factors:  (1)  The  energy  re¬ 
leased  by  the  explosion;  (2)  the  position  of  the  explosive  rela¬ 
tive  to- the  surface;  (3)  the  material  type;  and  (4)  gravitational 
effects.  The  influence  of  the  energy  release  is  obvious,  the 
larger  the  charge,  the  larger  the  crater.  When  the  charge  is 
on  or  above  the  ground  surface,  cratering  effects  are  small. 

As  the  charge  is  placed  deeper  in  the  ground,  the  size  of  the 
crater  increases,  both  in  radius  and  depth,  until  a  maximum  is 
reached  after  which  the  crater  size  will  decrease  with  in¬ 
creasing  depth  of  burial.  For  deeply  buried  explosives,  an 
underground  cavity  is  formed.  The  surface  itself  may  be  raised 
and  may  eventually  subside  to  form  a  depression  crater.  For 
materials  which  bulk  during  the  explosion  process,  there  is  a 
region  where  rubble  mounds  will  be  formed.  Mounding  and 
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subsidence  depend  upon  the  material  and  depth  of  burst.  The 
larger  the  material's  energy-dissipative  properties,  the 
smaller  the  crater  size  will  be.  The  size  of  the  apparent  crater 
is  affected  by  the  amount  of  fall-back  material,  that  is,  the 
material  originally  ejected  which  returns  under  gravity  to 
the  crater  zone. 

Crater  Terminology  (22,27).  A  few  basic  terms  and 
definitions  are  presented  to  provide  an  acquaintance  with 
significant  zones,  dimensions  and  terminology  used  throughout 
the  report.  Again  referring  to  Fig.  1,  the  cross  section  of  a 
typical  crater  and  the  adjacent  zones  of  disturbance  are  shown. 

The  apparent  crater  is  defined  as  that  portion  of  the 
visible  crater  which  is  below  the  preshot  ground  surface.  The 
apparent  crater  would  be  the  net  design  excavation  for  most 
engineering  applications. 

The  true  crater  is  defined  as  the  boundary  (below  preshot 
ground  level)  between  loose,  broken,  disarranged  fall-back 
materials  and  the  underlying  rupture  zone  material  which  has 
been  crushed  and  fractured,  but  has  not  experienced  significant 
displacement  or  disarrangement.  The  true  crater  boundary  is 
not  a  distinct  surface  of  discontinuity,  but  rather  a  zone  of 
transition  between  the  rupture  zone  and  fall -back  materials. 

The  apparent  lip  is  composed  of  two  parts,  the  true  lip, 
which  is  formed  by  the  upward  displacement  of  the  ground  surface, 
and  the  ejecta  material  deposited  on  the  true  lip. 
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The  visible  crater  comprises  the  apparent  crater  and  the 


apparent  lip. 

The  fall-back  consists  of  natural  materials  which  have 
experienced  significant  disarrangement  and  displacement  and 
have  come  to  rest  within  the  true  crater. 

The  ejecta  consists  of  material  thrown  out  above  and 
beyond  the  true  crater. 

The  rupture  zone  is  that  zone  extending  from  the  true 
crater  boundary  in  which  crushing  and  fracturing  have  occurred. 

The  plastic  zone  is  that  portion  of  the  cratered  medium 
beyond  the  rupture  zone  in  which  permanent  deformation  has 
occurred. 

The  elastic  zone  extends  beyond  the  plastic  zone  and  is 
characterized  by  the  absence  of  blast  produced  fissures, 
cracks,  or  permanent  displacement  of  material. 

The  optimum  depth  of  burst  is  that  depth  for  a  specified 
explosive  charge  which  produces  the  largest  crater.  For  any 
one  charge  in  any  one  medium,  there  may  be  three  optimum  depths 
of  burst  depending  upon  whether  the  largest  radius,  depth,  or 
volume  is  desired. 

A  scaled  dimension  refers  to  a  particular  crater  dimension 
divided  by  the  explosive  weight  of  the  charge  to  some  power 
(usually  between  1/4  and  1/3).  If  the  proper  exponent  is 
sele  ted,  scaling  between  different  sized  explosive  charges  for 
a  particular  scaled  depth  of  burst  is  then  possible. 
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Crater  dimensional  data  used  in  this  report  along  with 
accompanying  symbols  and  definitions  is  given  in  Fig.  2.  A 
detailed  description  of  all  pertinent  single-charge  crater 
dimensional  data  can  be  found  in  reference  22. 

Cratering  Mechanisms  (27,38,46).  The  primary  mechanisms 
or  processes  which  produce  craters  in  rock  or  soil  may  be 
categorized  as:  (1)  Vaporization  and  melting  of  the  material 
immediately  surrounding  the  source  of  a  nuclear  explosion; 

(2)  crushing,  compaction,  fracturing* and  plastic  deformation 
of  the  medium  closely  surrounding  the  explosive  gas  cavity; 

(3)  spalling  of  the  surface;  (4)  acceleration  of  the  fractured 
material  overlying  the  explosion  by  trapped  gases;  and  (5)  sub¬ 
sidence  and  fail-back  of  the  material  as  the  explosive  pressure 
goes  to  zero  and  the  force  of  gravity  predominates. 

The  tremendous  pressures  resulting  from  an  explosive 
detonation  (10-100  million  atmospheres  for  a  nuclear  explosive) 
generate  a  shock  wave  which  propagates  as  a  high-pressure 
discontinuity.  This  high-pressure  discontinuity,  or  shock  front, 
transfers  energy  to  the  medium,  and  in  turn,  alters  the  physical 
characteristics  of  the  medium.  In  the  immediate  vicinity  of  a 
nuclear  explosion,  vaporization  and  melting  of  the  material 
occurs.  The  peak  pressure  in  the  shock  wave  diverges  and  energy 
is  expended  in  doing  work  on  the  medium.  When  the  pressure 
and  shear  stress  levels  exceed  the  dynamic  crushing  strength  of 
the  material,  work  on  the  medium  is  manifested  in  crushing, 
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R  -  Radius  of  apparent  crater  measured  at  preshot 
ground  surface 

D  -  Maximum  depth  of  apparent  crater  below  preshot 
ground  surface 

H  i  -  Apparent  crater  lip  crest  height  above  preshot 
ground  surface 

Va  -  Volume  of  apparent  crater  below  preshot  ground 
surface 

Dob  -  Depth  of  burst  (distance  to  ZP  from  SGZ) 

ZP  -  Zero  Point-effective  center  of  explosion  energy 

SGZ  -  Surface  Ground  Zero  (point  on  surface  vertically 
above  ZP) 

FIG.  2.  DIMENSIONAL  DATA  FOR  SINGLE  CHARGE  CRATER 
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heating,  displacement, and  deformation  of  the  material.  When 
the  compressive  and  shear  waves  which  propagate  from  the 
detonation  encounter  the  surface,  a  tensile  wave  and  another 
shear  wave  are  reflected,  and  since  the  tensile  strengths  of 
rock  and  soil  are  much  less  than  their  compressive  strengths, 
spalling  of  the  surface  occurs.  Rock  also  tends  to  spall  along 
pre-existing  fractures  and  planes  of  weakness.  The  first  two 
processes  may  be  classified  as  short  term  mechanisms  since 
they  last  only  a  fraction  of  a  second.  Gas  acceleration,  on 
the  other  hand,  is  a  comparatively  long-period  process  which 
imparts  motion  to  the  material  around  the  explosion  by  the 
adiabatic  expansion  of  gases  trapped  in  the  cavity.  Finally 
the  force  of  gravity  pulls  all  the  overlying  fractured  and 
crushed  material  into  the  cavity  and  pulls  all  the  loose 
material  thrown  into  the  air  by  spalling  and  gas  acceleration 
back  into  and  around  the  crater.  What  remains  is  an  apparent 
crater  underlayed  with  crushed  and  displaced  material. 

Effects  of  Depth  of  Burst.  The  part  each  of  the  above 
mechanisms  play  in  producing  a  crater  is  very  strongly 
dependent  upon  the  scaled  depth  of  burst  of  the  explosion  and 
the  medium  in  which  the  detonation  occurs.  Shown  in  Fig,  ? 
are  typical  crater  cross  sections  in  rock  showing  the  effect 
of  depth  of  burst.  As  summarized  by  Nordyke  (46),  Fig.  4, 
the  contribution  of  each  of  the  mechanisms  to  apparent 
crater  depth  as  a  function  of  the  charge  depth  of  burst  js  shown. 
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SECTION  III 


REVIEW  OF  PREVIOUS  CRATERING  RESEARCH 

There  have  been  thousands  of  cratering  experiments  and 
hundreds  of  technical  reports,  articles, and  books  written  as  a 
result  of  these  experiments  and  associated  phenomena.  Over  300 
of  these  publications  were  reviewed  for  this  study.  This 
section  presents  a  brief  review  of  cratering  experience  and 
analysis.  A  brief  history  of  recent  cratering  research  and 
symposiums  relating  to  that  research  is  presented  followed  by 
the  various  prediction  techniques  available  and  the  information 
found  on  the  effects  of  material  properties. 

Nuclear  Cratering  Experience  (13,45).  There  have  been  10 
nuclear  detonations  at  the  Nevada  Test  Site  involving  four 
different  geologic  media  in  level  terrain  which  resulted  in  the 
creation  of  craters  (or  a  mound  as  in  the  case  of  Sulky)  which 
are  considered  applicable  to  explosive  excavation.  These  de¬ 
tonations  are  listed  as  the  first  10  events  in  Appendix  I. 
Yields  for  these  events  varied  from  0.42  kiloton  for  Danny  Boy 
to  100  kilotons  for  Sedan  (on o  kiloton  =  1,000  tons  equivalent 
weight  of  TNT). 

High-Explosive  Cratering  Experience.  Good  summaries  of 
past  research  in  this  area  can  be  found  in  references  46,  81 
and  87.  Lists  for  the  majority  of  the  experiments  can  be  found 
in  references  15,  60  and  87.  In  all,  these  experiments  number 
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in  the  thousands  and  involve  more  than  20  different  soil  and  rock 
materials.  They  range  in  yield  from  1  gram  to  1  million  pounds 
of  TNT  and  include  the  full  range  of  depths  of  burst.  Con¬ 
sidering  each  material,  the  largest  number  of  these  experiments 
were  performed  in  a  lightly  cemented  sand-gravel  mixture  known 
as  desert  alluvium.  The  number  of  experiments  applicable 
for  this  research  was  considerably  less  than  the  total  number 
available.  For  instance,  original  plans  had  included  the 
Dugway  cratering  series  in  limestone,  granite, and  sandstone  (76). 
Later  review,  however,  determined  that  the  crater  dimensions 
reported  were  not  those  for  the  apparent  crater  as  had  been 
indicated  by  Vortman  (87),  but  were  mucked  dimensions;  i.e., 
all  loose  rubble  in  the  crater  was  removed  before  crater  dimen¬ 
sions  were  measured.  Crater  dimensions  for  this  series,  there¬ 
fore,  lie  somewhere  between  apparent  and  true  crater  definitions. 
It  was  also  necessary  to  eliminate  hundreds  of  experiments 
reported  by  Sager  (60)  because  other  than  spherical  charges  were 
used  and  because  very  few  of  the  experimentors  had  reported 
material  properties.  A  good  conclusion  that  can  be  made  regarding 
the  various  experiments  and  their  data  plots  is  that  the  data 
are  highly  scattered  because  explosive  experiments  in  natural 
materials  are  difficult  to  control.  The  high-explosive  experi¬ 
ments  used  in  this  research  are  shown  in  Appendix  I  beginning 
with  event  9. 

Symposiums  Related  to  Cratering.  The  Plowshare  Program 
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(the  peaceful  uses  of  nuclear  explosives)  was  formally  estab¬ 
lished  in  1957  (30).  At  about  the  same  time,  the  first  con¬ 
tained  nuclear  experiment,  Ranier,was  executed  (31).  Ranier 
was  a  1.7-kiloton  explosion  900  feet  below  the  surface.  Be¬ 
cause  of  the  success  of  the  Ranier  event,  numerous  speculations 
were  made  as  to  the  uses  of  underground  explosions.  The  gen¬ 
eral  range  of  ideas  was  first  reported  publicly  in  1958  at 
the  Atoms  for  Peace  Conference  in  Geneva.  This  later  became 
the  First  Plowshare  Symposium.  In  1959,  the  Second  Plowshare 
Symposium  was  held  in  San  Francisco  and  initial  presentations 
on  the  uses  of  nuclear  explosives  for  excavation  and  extrapo¬ 
lation  of  chemical  explosive  data  for  possible  nuclear  explosive 
application  were  presented  (54). 

In  1961,  a  Geophysical  Laboratory  -  Lawrence  Radiation 
Cratering  Symposium  was  held  in  Washington  D.C.  (44).  This  was 
really  the  first  somewhat  all  inclusive  up-to-date  presentation 
relative  to  cratering  data  and  phenomena.  Scaling  laws, 
empirical  analysis,  theoretical  calculations,  nuclear  cratering 
to-date, and  explosive  craters  in  desert  alluvium,  tuff  and 
basalt  were  some  of  the  more  pertinent  subjects  presented. 
Nordyke  presented  his  preliminary  theory  for  the  mechanics  of 
crater  formation;  this  theory  is  still  being  followed  today. 

In  1964,  the  Third  Plowshare  Symposium  with  its  theme 
"Engineering  with  Nuclear  Explosives"  was  held  at  the  University 
of  California,  Davis  Campus  (74).  Over  30  presentations  were 
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made  that  offered  an  up-to-date  picture  for  using  nuclear 
explosives  for  engineering  purposes.  By  this  time  additional 
nuclear  cratering  data,  including  the  100-kiloton  Sedan  event, 
were  available.  Knox  and  Terhune  (32)  presented  results  ob¬ 
tained  *rom  using  SOC,  the  two-dimensional  computer  model  of 
cratering  physics  during  the  gas  acceleration  phase.  A  good 
portion  of  the  data  presented  at  the  conference  was  later  used 
by  Teller  (70)  to  write  the  book  The  Constructive  Uses  of 
Nuclear  Weapons. 

In  January  1970,  the  most  recent  symposium,  Engineering 
with  Nuclear  Explosives,  was  held  in  Las  Vegas,  Nevada  (55). 

Over  100  presentations  were  made  of  which  17  were  directly 
related  to  excavation.  Since  a  number  of  the  presentations  re¬ 
late  to  this  research,  they  are  discussed  in  the  next  two  sub¬ 
sections  . 

Prediction  Techniques  and  Scaling.  To  date  there  are 
basically  two  approaches  being  used  for  predicting  crater 
dimensions.  The  first  involves  computer  calculations  of  the 
mound  and  cavity  growth  used  in  conjunction  with  a  free-fall, 
throw  out  model  which  gives  an  estimate  of  the  crater  radius  and 
ejecta  boundary.  The  second  basic  approach  uses  empirical 
scaling  relationships.  Another  approach,  although  not  so  widely 
used  but  which  deserves  mentioning,  is  a  quasi -static  approach 
which  considers  cratering  as  an  earth  pressure  problem. 

SOC  (spherical  symmetry,  one  dimensional)  and  TENSOR 
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(cylindrical  symmetry,  two  dimensional)  computer  codes  numeri¬ 
cally  describe  the  propagation  of  a  stress  wave  of  arbitrary 
amplitude  through  a  medium  (6,12,32,36,37,71).  These  codes  are 
hydrodynamic  Lagrangian  finite-difference  approximations  of  the 
equation  of  motion  which  describe  the  behavior  of  a  medium 
subjected  to  a  stress  tensor  in  one  (SOC)  ana  two  (TENSOR) 
space  dimensions.  The  code  calculations  handle  both  the  initial 
shock  wave  phase,  which  creates  spall  velocities  and  the  gas 
acceleration  phase.  The  end  product  of  the  TENSOR  code  cal¬ 
culations  is  a  chronological  history  of  the  cavity  and  mound 
growth  resulting  from  an  underground  explosive  detonation.  The 
code  calculations  runs  until  the  particle  velocities  no  longer 
increase  significantly  from  cycle  to  cycle.  At^this  point,  a 

free-fall  throw  out  model  calculation  is  used  to  ^determine  the 

\ 

mode  of  deposition  of  that  material  which  has  beenv- given  suf¬ 
ficient  velocity  to  pass  the  original  ground  surface.  The 
ballistic  trajectory  of  any  given  mass  determines  its  final 
position  on  the  surface.  The  throw  out  model  calculation 
permits  one  to  estimate  crater  radius  and  the  maximum  range  to 
which  significant  material  is  thrown  by  the  detonation.  An 
estimate  of  the  crater  depth  may  also  be  made  by  considering 
the  stability  of  the  cavity  walls  and  the  bulking  characteristics 
of  the  material  which  falls  back  into  the  crater  opening.  Since 
these  codes  assume  that  the  material  behaves  hydrodynamically 
not  only  in  the  melted  region  but  in  the  consolidated  and 
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cracked  regions  of  the  media,  extensive  laboratory  testing  is 
required  to  develop  pressure-volume  relationships  for  the 
material  involved.  The  big  advantage  to  this  method  is  that 
it  provides  a  visual  graphical  representation  of  the  mound 
and  cavity  growth  and  fall-back.  The  disadvantages  are  that 
it  requires  an  enormous  amount  of  computer  time  to  accomplish 
the  computation  and  it  requires  extensive  material  testing  to 
obtain  pressure-volume  relationships,  rigidity  modulus,  ten¬ 
sile  strength  and  distortional  energy  limits.  Failure  in  the 
material  is  assumed  to  occur  either  when  the  tensile  strength 
is  exceeded  or  when  maximum  distortional  strain  energy  is 
reached. 

The  second  crater  geometry  prediction  approach  involves 
the  use  of  scaling  laws  which  relate  crater  dimensions  from 
some  reference  yield  to  crater  dimensions  for  any  energy  yield. 
The  preponderance  of  cratering  research  (8,10,11 ,15,46,61 ,77,88) 
has  been  concerned  with  or  utilized  this  prediction  approach. 
This  procedure  in  its  simplest  form  is  based  ori  an  empirical 
determination  of  the  scaling  exponent,  m,  as  a  function  of  soil 
type,  using  the  assumed  relationships 

R  =  CnW1/m;  D  =  C2W1/m;  V  =  C3W1/m . (1) 

where  R,  0  and  W  are  the  radius,  depth  and  volume,  respecti <"=’y, 
of  the  crater,  are  constants  related  to  the  soil  type  and 
depth  of  burst  and  W  is  the  energy  release.  If  a  dimensional 
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analysis  is  made,  it  would  appear  that  m  should  be  3  for  radius 
and  depth  and  1  for  the  volume  if  the  effects  of  gravity  and 
material  friction  are  omitted;  however  when  gravity  of  the 
material  is  considered,  m  becomes  4  for  radius  and  depth  and  4/3 
for  volume.  The  results  of  cratering  experiments  to  date, 
however,  have  led  to  the  development  of  m  =  3.4  for  radius  and 
depth  and  m  =  3.4/3  for  volume.  Although  this  relationship  is 
very  simple,  the  only  way  can  be  determined  is  by  performing 
a  series  of  cratering  experiments  for  the  material  in  question. 
There  is  considerable  actual  data  scatter,  however,  compared 
with  the  smooth  curve  this  technique  produces.  The  constants, 

C.j ,  therefore , appear  to  be  variables  which  depends  on  material 
properties.  At  best,  then,  these  relationships  provide  only  a 
very  rough  estimate  of  cratering  dimensions  since  all  material 
properties  are  ignored. 

A  third  method  for  predicting  crater  dimensions  considers 
cratering  as  an  earth  pressure  problem  (42,78,79,80,81).  This 
method  considers  expansion  of  the  explosive  cavity  to  some 
ultimate  radius  and  pressure,  at  which  time  equilibrium  is  as¬ 
sumed  to  exist,  at  least  for  a  while.  For  the  cratered  medium, 
it  is  assumed  that  a  part  of  it,  adjacent  to  the  cavity, 
behaves  as  a  rigid  plastic  solid,  defined  by  a  Mohr's  envelop  for 
the  material.  At  a  sufficient  distance  from  the  explosive 
charge  the  medium  is  assumed  to  behave  as  a  linearly  deformable, 
isotropic  solid  defined  by  a  deformation  modulus  and  a  Poisson's 
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ratio.  This  method  appears  very  reasonable,  however  it  has 
only  been  applied  to  very  small  scale  laboratory  experiments. 

In  addition,  it  does  not  consider  surface  spallation  and  does 
not  seem  to  apply  for  surface  bursts. 

Material  Property  Effects.  Although  numerous  observations 
and  studies  of  cratering  for  both  conventional  and  nuclear 
explosives  have  led  to  a  fair  understanding  of  the  basic 
processes  and  phenomena  involved, the  physical  characteristics 
of  the  cratering  medium  which  significantly  influence  crater 
sizes  and  shapes  were  found  to  be  largely  unknown. 

Whitman  (89)  postulated  that  crater  dimensions  were 
primarily  related  to  soil  type  through  the  soil's  shearing 
resistance.  His  study  developed  trends  between  eraser  size 
and  soil  strength  for  the  weaker  soils,  but  trends  for  the 
stronger  soils  and  rocks  were  quite  obscure. 

Based  on  limited  data,  Baker  (1)  was  able  to  relate  material 
seismic  velocity,  angle  of  shearing  resistance  for  sands,  tensile 
strength  for  rocks, and  relative  consistency  of  clay  to  a  radius 
modulus  developed  by  Saxe  (6?)  to  relate  a  scaled  normalized 
radius  and  the  scaled  depth  o*  burial. 

Chabai  (11)  made  a  dimensional  analysis  of  scaling 
dimensions  of  craters  and  considered  the  following  medium 
properties  as  being  sufficient  to  describe  the  phenomena  of 
cratering:  density  of  undisturbed  medium,  yield  strength  of 
the  medium,  a  viscosity  or  dissipation  variable  of  the  medium 
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and  the  sonic  velocity  in  the  medium. 

Westine  (88)  considered  dimensional  analysis  and  developed 
the  following  relationship: 

Ra/Dob  =  f[W 7/24/(p7/24  g1/8  c1/3  Dob)]  ....  (2) 
in  which  p  =  mass  unit  weight,  g  =  acceleration  due  to  gravity 
and  c  =  seismic  velocity.  Wes  tine's  plots  of  the  data  from 
five  different  materials  show  some  reduction  of  the  data 
scatter  using  his  technique.  The  disadvantage  in  using  this 
approach  is  that  it  implicitly  assumes  that  an  increase  in  unit 
weight  and  seismic  velocity  will  result  in  a  smaller  crater 
dimension. 

Terhune  (71)  listed  the  following  material  parameters  in 
the  order  of  their  importance  for  determining  the  cratering 
efficiency  of  the  medium:  (1)  Water  content;  (2)  shear 
strength;  (3)  porosity  (compactibility);  and  (4)  compressibility. 
He  stressed  the  point  that  an  increase  in  water  content  decreases 
the  compressibility,  drastically  reduces  the  shear  strength  and 
provides  an  additional  energy  source  in  the  form  of  a  non¬ 
condensable  gas. 
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SECTION  IV 


PURPOSE,  SCOPE, AND  PROCEDURE 

The  purpose  of  this  research  was  to  evaluate  the  data  from 
published  cratering  experiments  in  an  effort  to  show  the  effect 
of  soil  and  rock  properties  on  crater  dimensions  in  conjunction 
with  the  burial  depth  and  energy  of  the  explosive  charge.  A 
secondary  purpose  was  to  show  that  no  real  analysis  can  ever 
be  made  of  cratering  data  until  controlled  experiments  are  con¬ 
ducted  or  unless  soil  and  rock  property  measurements  are  care¬ 
fully  made  throughout  the  material  field. 

Of  the  thousands  of  cratering  experiments  that  have  been 
conducted,  only  a  little  more  than  200  tests  were  selected  for 
analysis.  The  remaining  tests  did  not  meet  the  general  criteria 
established  for  this  study.  From  the  tests  selected,  those  where 
the  charge  was  detonated  above  the  surface  or  the  depth  of  burial 
was  so  great  that  a  mound  or  slight  depression  developed,  instead 
of  a  crater,  were  not  used  in  the  analysis. 

The  procedure  used  to  accomplish  the  research  was 
statistical  analysis  of  data  from  existing  cratering  experiments. 
This  involved  the  following:  cataloging  all  crater  and  related 
material  property  data;  selecting  a  proper  regression  model  to 
develop  empirical  equations  for  apparent  crater  radius,  depth 
and  volume  in  terms  of  explosive  weight  and  depth  of  burst  for 
specific  materials;  integration  of  the  various  material  properties 
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into  the  model  to  develop  the  best  equations  for  all  the 
materials  combined;  and  analyzing  the  best  equations  to 
determine  the  specific  effect  of  the  material  properties  used. 
An  inherent  part  of  the  procedure  was  the  development  and  use 
of  computer  programs. 


SECTION  V 


EXPERIMENTAL  DATA 

This  section  presents  the  method  and  approach  used  to 
acquire  and  catalog  the  experimental  crater  and  associated 
material  property  data  necessary  for  this  study.  The  cataloged 
listing  is  referred  to  and  discussed  in  some  detail.  Sorting 
of  the  data  to  allow  for  better  comparison  and  analysis  is 
briefly  described. 

Data  Acquisition.  Using  the  Corps  of  Engineers'  "Com¬ 
pendium  of  Crater  Data"  (60),  Vortman's  "Ten  Years  of  High 
Explosive  Cratering  Research  at  Sandia  Laboratory"  (8/),  and 
Circeo's  "Nuclear  Excavation:  Review  and  Analysis"  (13),  as 
guides  to  previously  conducted  experiments,  the  original  source 
documents  were  obtained  and  reviewed.  This  involved  making 
visits  to  Sandia  Laboratory  and  Air  Force  Weapons  Lab,  Albuquerque, 
N.M.  and  to  Lawrence  Radiation  Laboratory  and  U.  S.  Army  Engineer 
Nuclear  Cratering  Group,  Livermore,  California. 

A  review  of  the  literature  showed  that  there  was  a  wide 
variation  in  the  experiments  as  well  as  a  wide  variation  in  the 
parameters  measured.  The  following  criteria  was  established  to 
determine  if  a  particular  cratering  experiment  was  to  be 
cataloged  and  used  in  this  analysis: 

1.  The  explosive  charge  had  to  be  single  and  spherical 
with  a  TNT  equivalent  weight  of  at  least  1  pound. 
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2.  The  dimensions  measured  had  to  be  for  the  apparent 
crater. 

3.  The  experimental  terrain  had  to  be  reasonably  level. 

4.  Sufficient  material  properties  had  to  be  measured  or 
it  had  to  be  possible  to  estimate  them  with  some 
degree  of  confidence  from  other  recorded  data. 

Although  crater  dimensions  were  obtained  for  the  various 
experiments  with  relative  ease,  soil  and  rock  properties 
presented  a  very  perplexing  problem.  For  a  number  of  the 
nuclear  events,  extensive  soil  borings  and  tests  were  reported. 
In  other  cases,  only  one  soil  test  pit  and  limited  testing  was 
reported  for  a  complete  series  of  experiments.  Somewhat 
arbitrarily,  but  keeping  in  mind  the  cratering  phenomenology 
involved,  those  material  properties  which  were  thought  to 
effect  crater  geometry  were  selected  to  be  cataloged. 

Data  Cataloging.  By  considering  all  the  data  required  to 
properly  catalog  each  cratering  event,  formats  and  programs  for 
computer  input  and  output  were  developed.  Every  effort  was 
made  to  keep  the  data  for  a  particular  cratering  event  to  a 
minimum  but  yet  make  the  data  as  complete  as  possible.  For 
example,  cataloging  only  the  grain  specific  gravity,  dry  unit 
weight  and  moisture  content  of  the  media  for  a  particular  event 
allowed  for  computer  calculation  of  various  parameters  such  as 
total  unit  weight,  degree  of  saturation,  porosity,  percent 
air,  void  ratio,  etc.  Data  input  to  the  computer  was  made  to 
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serve  a  dual  purpose.  It  was  used  to  produce  the  cataloged 
listing  of  crater  data  and  as  a  basic  data  source  for  the 
analysis. 

The  listing  of  all  data  cataloged  for  this  research  is 
included  in  Appendix  I.  This  appendix  is  divided  into  the 
following  four  segments: 

1.  A  list  of  all  notation  and  definitions  associated 
with  the  cataloged  data. 

2.  The  computer  output  crater  data  list. 

3.  The  two  line  computer  listing  of  all  the  cataloged 

material  property  data  associated  with  the  crater 
data  list. 

4.  The  notes  referred  to  in  both  the  crater  data  and 
material  property  data  listings. 

It  was  initially  thought  that  there  existed  many  more 
measurements  of  cratering  event  material  property  data.  In 
the  final  analysis,  sufficient  data  was  just  not  avail¬ 
able  to  obtain  measured  values  for  each  and  every  event  or 
in  most  cases  even  for  a  series  of  events.  Estimated  values 
dominate  the  material  property  data  section  of  the  cataloged 
data.  To  differentiate,  estimated  values  are  indicated  by  a 
dollar  sign  in  the  data  listing. 

Viewing  the  crater  data  list,  it  can  be  seen  that  the 
approximate  energy  of  the  nuclear  explosive  cratering  events 
was  only  estimated  within  20  percent.  This  was  due,  apparently. 
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to  the  unsureness  of  the  amount  of  nuclear  material  which 
actually  participated  in  the  reaction.  As  discussed  by 
Vortman  (87)  for  most  high  explosive  events,  cast  spherical 
charges  of  TNT  detonated  at  their  centers  were  used  for  yields 
of  1000  pounds  or  less.  For  charges  larger  than  1000  pounds, 
cast  blocks  of  TNT  were  stacked  to  approximate  a  sphere. 

For  certain  experiments  liquid  nitromethane  was  used.  The 
liquid  was  placed  in  a  spun  aluminum  sphere  or  in  a  mined 
cavity  lined  with  an  impervious  material.  For  a  number  of 
the  very  small  explosive  tests,  Military  C-4  explosive  was 
used.  Equivalent  weight  of  TNT  factors  for  these  latter 
explosives  was  based  on  the  heat  of  detonation  as  reflected 
in  Cook  (14).  There  is,  however,  some  question  as  to  whether 
these  equivalency  factors  are  completely  valid.  There  seems 
to  be  an  optimum  rate  of  burning  for  an  explosive  for  a 
particular  depth  of  burst  in  a  particular  material  which  will 
produce  the  largest  crater.  It  is  felt,  however,  that  thfj 
equivalent  weight  of  TNT  shown  in  the  crater  data  list  for  all 
events  except  the  nuclear  tests  is  within  five  percent. 

Depending  on  the  method  used  by  the  original  investigators 
to  measure  crater  dimensions,  these  dimensions  are  considered 
accurate  only  to  within  five  percent.  Again  as  discussed  by 
Vortman  (87),  crater  measurements  have  been  determined  using 
various  techniques.  These  techniques  have  consisted  of  almost 
everything  imaginable  from  conventional  ground  surveys  to  the 


26 


use  of  adjustable  rods  to  aerial  mapping.  Lip  heights  and 
slope  angles  shown  in  the  listing  (although  not  specifically 
used  in  this  research)  were  not  reported  by  many  of  the 
investigators.  Values  reported,  in  many  cases,  reflect  values 
measured  from  typical  cross  sections. 

Material  property  values  reported  in  the  material  property 
data  list  reflect  those  which  were  either  reported  by  the 
investigator  as  an  average  for  the  experiment  or  series  of 
experiments,  or  where  possible,  taken  from  boring  logs  and 
associated  tests.  In  the  latter  case,  a  weighted  average  for 
a  particular  material  property  was  computed  for  the  vertical 
column  in  question.  These  weighted  averages  for  the  various 
borings  were  then  interpolated  or  extrapolated  to  obtain  the 
values  reported  for  the  explosive  event  location.  Although 
the  material  property  data  list  contains  predominately 
estimated  values,  the  majority  of  these  estimated  values  were 
extrapolated  from  experiments  in  like  or  similar  materials. 

It  would  have  been  advantageous  to  have  been  able  to  obtain 
specific  measured  values  for  each  event.  Where  measured, 
material  property  values  are  considered  accurate  only  to 
within  10  percent.  Where  values  are  estimated,  it  is  believed 
that  they  are  within  20  percent. 

Data  Sorting.  To  allow  for  better  comparison  and  analysis 
of  the  data  cataloged,  a  computer  subroutine  was  written  and 
used.  This  subroutine  provided  for  sorting  on  any  three 
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specified  parameters  at  one  time.  The  sorting  scheme  that 
was  found  to  be  most  useful  aligned  the  data  by  type  of 
material,  then  by  scaled  depth  of  burst  and  lastly  by  explosive 
weight. 
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SECTION  VI 


EMPIRICAL  DATA  ANALYSIS 

This  section  describes  the  regression  analysis  technique 
used  to  develop  a  functional  relationship  between  crater  dimen¬ 
sions  and  explosive  charge  weight  and  depth  of  burst  for  a  given 
soil  or  rock  type.  It  also  discusses  the  function  constants 
obtained,  presents  plots  of  the  actual  data  and  presents  the 
resultant  regression  curves  to  predict  crater  dimensions  for 
nine  groupings  of  the  data. 

Regression  Analysis  (16,39).  The  data  for  this  research 
were  analyzed  and  the  prediction  formulas  for  crater  dimensions 
were  developed  using  applied  regression  analysis.  First,  the 
dependent  variable  (the  variable  for  which  prediction  was 
desired)  was  selected  and  the  independent  variables  (the 
parameters  thought  to  have  some  influence  on  the  dependent 
variable)  were  assumed.  Next  a  form  of  the  answer  (the  model) 
was  assumed  and  the  data  applied  to  the  model  to  obtain  its 
coefficients.  This  was  accomplished  through  the  method  of 
least  squares  surface  fitting,  whereby  the  sum  of  squares  of 
the  distances  between  the  assumed  surface  and  the  actual  data 
points  was  minimized.  If  the  sum  of  squares  is  a  minimum  and 
the  coefficients  of  the  assumed  model  are  determined,  then  this 
is  the  best  fitting  surface  for  the  model  and  data  used.  Lastly 
through  computation  of  a  multiple  correlation  coefficient  and 
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by  an  examination  of  the  least  squares  residuals  the  prediction 
formula  was  evaluated. 

Regression  analysis  has  become  a  fairly  common  tool  for 
analyzing  experimental  data  and  developing  function  relation¬ 
ships  for  that  data.  For  this  reason,  it  seems  sufficient  to 
mention  only  that  an  enormous  amount  of  mathematical  computa¬ 
tion,  including  the  solution  of  a  large  number  of  simultaneous 
equations,  is  required.  Only  through  the  use  of  a  high-speed 
computer  is  this  possible. 

The  computer  program  used  was  specifically  written  for 
this  research.  Although  library  regression  analysis  programs 
were  available,  it  was  felt  that  these  programs  were  too  elabo¬ 
rate  and  did  not  provide  the  flexibility  needed  for  this 
research.  The  computer  program  written  was  intended  to  be  a  very 
flexible,  minimum  essential  program  that  would  adequately 
accomplish  the  data  cataloging  and  sorting,  the  least  squares 
surface  fit  and  the  evaluation  of  the  fit.  The  multiple  linear 
regression  portions  of  the  program  were  written  using  statistics 
books  as  guides  (16,24).  The  matrix  inversion,  multiplication 
and  print  subroutines  were  available  from  previous  research  (75) 
and  were  modified  for  use  here.  As  the  research  progressed, 
numerous  changes,  additions  and  deletions  were  made  based  on 
the  scheme,  procedure  or  purpose  being  tried  at  the  moment. 

The  program  was  originally  written  for  the  IBM  360/65  computer, 
but  was  later  revised  for  use  on  the  CDC  6600.  Appendix  II 
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contains  a  brief  description  of  the  program  and  its  essential 
features  along  with  a  typical  print-out  of  the  essential 
portions  of  the  program  and  its  output  after  being  run  on  the 
CDC  6600. 

Regression  Approaches  Considered  and  Tried.  The  number 
of  regression  approaches  considered  and  tried  was  so  numerous, 
it  is  superfluous  to  list  them  all.  At  the  beginning,  regres¬ 
sions  of  radius,  depth,  or  volume  as  a  linear  function  of  depth 
of  burst,  explosive  weight  and  various  material  properties  were 
attempted.  A  second  approach  considered  a  linear  crater  dimen¬ 
sion  to  be  a  function  of  various  dimensionless  parameters  taken 
from  Chabai *s  work  (11).  An  attempt  was  also  made  to  consider 
all  material  parameters  which  would  reflect  the  amount  of  energy 
being  dissipated  during  the  cratering  process.  A  closely 
related  approach  was  to  consider  those  material  properties  which 
would  relate  to  the  mechanisms  of  compaction,  subsidence,  spall 
and  gas  acceleration  as  proposed  by  Nordyke  (46).  In  all  of  the 
above  cases,  sufficient  material  properties  were  not  available 
to  include  all  applicable  terms  thought  to  be  important.  An 
attempt  was  made  to  use  available  and  applicable  material 
parameters  in  a  linear  fit,  but  to  no  avail.  When  one  considered 
a  second  order  regression  model,  the  number  of  parameters  to  be 
used  Suddenly  becomes  excessive  for  the  data  available  and  for 
the  basic  research  purpose. 

Of  special  note  wa s  the  attempt  to  use  Westine's  technique 
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(88).  Regression  using  this  approach  produced  a  high  multiple 
correlation  ratio.  However,  the  residuals  between  the  estimated 
and  actual  values  of  the  dependent  variables  were  found  to  be 
excessive.  This  was  particularly  true  for  the  range  of  data 
where  R  /Dob  was  loss  than  2. 

a 

Regression  Model  Used.  After  much  deliberation  and  due 
consideration  of  the  literature  reviewed,  it  was  felt  that  the 
better  approach  for  the  solution  of  crater  dimensions  for  one 
particular  media  was  to  consider  a  scaled  linear  crater  dimen¬ 
sion  as  a  function  of  the  scaled  depth  of  burst.  Although  there 
seemed  to  be  some  question  as  to  the  proper  scaling  to  use,  the 
scaling  exponent  which  appeared  to  most  recently  be  used  and 
justified  was  7/24.  This  figure  is  the  average  of  conventional 
cube  root  scaling  (gravity  effects  excluded)  and  fourth  root 
scaling  (gravity  effects  included).  The  scaling  exponent 
eventually  used  was  5/16  and  resulted  from  a  special  study  of 
the  data  cataloged  for  this  research.  The  next  question  which 
arose  was  what  model  should  be  used  for  regression?  After 
studying  plots  of  the  scaled  data  and  after  considering  what 
happens  to  crater  size  as  the  depth  of  burst  (or  height  of  burst) 
is  varied  from  one  extreme  to  the  other  it  became  obvious  that  a 
curve  reflecting  the  final  dimension  of  a  crater  when  plotted 
as  a  function  of  depth  of  burst  (and  height  of  burst)  should  be 
asymptotic  at  both  extremes  and  reach  its  maximum  value  at  the 
optimum  depth  of  burst.  This  suggested  the  bell  shaped  curve 
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with  its  inherent  advantages  and  disadvantages  which  will  be 
discussed  later.  Using  this  approach,  the  regression  suddenly 
improved  for  any  one  particular  cratering  media. 

The  general  form  of  a  bell  curve  is  as  follows: 

y  =  B1  exp  [B2  (x+B3)2] . (3) 

That  of  the  skewed  bell  curve  used  took  the  following  form: 

y  =  B-j  exp  [B2  (x+B3) (x+B4)2] . (4) 

in  which  x  =  the  scaled  depth  of  burst  (Dob/W^^)  and  where  y  = 
the  scaled  linear  apparent  crater  dimension  being  considered: 

(1)  Radius  (R^/W5/16) ;  (2)  depth  (D  /W5/16);  or  (3)  cube  root 
of  volume  (V^/W^^).  A  big  advantage  to  these  curves  is  that 
B^  is  the  maximum  height  of  the  curve  and  it  occurs  along  the 
abscissa  at  -B3  for  the  standard  curve  and  at  -B^  for  the 
skewed  form.  In  other  words  if  y  is  equal  to  the  scaled  radius 
for  the  first  case,  then  the  maximum  scaled  radius  is  B-j ,  and 
it  occurs  at  an  optimum  scaled  depth  of  burst  of  -B3>  sets 
the  rate  of  change  of  the  slope  away  from  the  (B^B^)  point. 

The  (x+B3)  term  in  the  second  case  allows  for  skewing  the 
standard  bell  curve.  As  can  be  seen,  this  introduces  a  second 
root  to  the  equation.  In  tne  majority  of  surface  fitting  cases, 
this  posed  no  problem  since  the  second  root  and  associated  slope 
change  were  outside  the  range  of  the  data.  Wh^re  the  data  was 
minimal  and  somewhat  scattered,  using  the  skewed  bell  curve 
proved  infeasible. 
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Since  the  above  curves  are  not  applicable  to  multiple 
linear  regression,  they  were  used  in  the  following  coded  forms. 
For  the  standard  bell  curve: 

In  y  =  A.|  +  A2x+  A3  x2  . (5) 

2 

where  A^  =  In  B-|  +  B2  B3  ,  A2  -  2B„  B^  and  A^  =  B2;  or  con¬ 
versely:  B-j  =  exp  (A-j  -  A22/4A3),  B2  =  A3  and  B3  =  A2/2A3> 

For  the  skewed  bell  curve: 

In  y  =  A-j  +  A2  x  +  A3  x2  +  A4  x2  .  .  .  (6) 

where  A^  =  In  B^  +  B2  B3  A2  =  B2  B^  (2B3+B^),  A3  =  B2 
(B3+2B^)  and  A^  =  B^ ;  or  conversely  B2  =  A^ ,  B^  =  [A0/A^  ± 
(A32/A42-3A2/A4)1/2]/3,  R?  =  A3/A,j  -  2B4  and  B]  =  exp(ArB2 
b3  B42) 

As  can  be  seen  for  this  later  case,  two  sets  of  values  for  the 
constants  B4,  B3  and  B^  are  obtained.  Examination  of  the 
numerical  values  obtained  for  a  particular  curve  fit,  however, 
quickly  show  the  correct  set  to  be  used. 

Determination  of  Best  Scaling  Exponent.  r o,  all  the 
initial  surface  fits  made,  1/4,  7/24  and  1/3  were  used  as  the 
scaling  exporents.  For  every  set  of  data  except  one,  7/24  was 
found  to  produce  the  best  fit  of  the  data.  Where  the  data 
consisted  of  a  large  amount  of  surface  bursts,  1/3  was  found 
to  produce  the  best  fit.  This  suggested  using  a  variable 
exponent  as  a  function  of  depth  of  burst.  Numerous  computer 
runs  v/ere  made  attempting  to  use  an  exponent  which  was  a 
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function  of  depth  of  burst  or  an  initial  scaled  depth  of  burst.  j 

Although  a  better  fit  was  obtained  for  surface  burst  data,  the  | 

over-all  fit  was  never  as  good  as  that  obtained  from  just  using  J 

7/24.  Vortman  (84)  showed  that  the  scaling  exponent  for  surface  j 

1 

bursts  could  vary  from  low  of  0.23  for  depth  to  0.44  for  radius 
depending  on  the  material  in  question. 

A  special  computer  run  was  made  to  vary  the  scaling 
exponent  from  a  value  of  0.29  to  0.35.  It  was  found  that  a 
scaling  exponent  of  0.31  produced  the  best  fits  for  radius  and  j 

cube  root  of  volume  while  0.32  produced  the  best  fit  for  depth 
(there  was  very  little  difference,  however,  between  the  fits 
obtained  from  using  0.31  and  0.32).  It  was  decided  to  use  : 

5/16  (or  0.3125)  as  the  scaling  exponent  for  the  remainder  of 
the  research.  This  figure  represents  the  average  of  the  7/24 
and  1/3  figures  which  have  been  used  so  predominately  by  other  > 

investigators.  , 

Equations  and  Daca  Plots  for  Specific  Materials.  Using 
equations  (3)  and  (4)  as  the  regression  models  (the  standard  •! 

i 

bell  curve  and  the  skewed  bell  curve),  the  coefficients  (B  ) 

( 

values)  were  generated  to  obtain  the  best  fit  of  the  data  to  \ 

^mDirically  predict  the  scaled  depth,  radius  or  cube  root  of 

volume.  This  was  accomplished  for  each  of  eight  groups  of 

material  (or  experiment)  data  and  in  addition  for  all  the  data 

combined.  Linear  dimensions  used  in  these  equations  were  in  feet 

while  W  was  in  pounds  of  TNT.  Listed  in  Table  1  are  the  B  values 
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EQUATION  COEFFICIENTS  FOR  VARIOUS  MATERIALS 


obtained.  The  number  of  8  values  listed  for  a  particular  crater 
dimension  for  a  particular  material  reflects  which  equation 
produced  the  best  fit.  Actual  plots  of  the  data  used  and  the 
resultant  empirical  curves  obtained  are  included  as  Appendix 
III. 

A  study  of  these  data  plots  reveals  several  interesting 
facts.  The  smallest  craters  were  produced  in  playa,  the 
lighter  weight,  weaker  material;  next  came  the  rocks,  alluvium 
and  sand;  and  lastly  the  largest  craters  were  produced  in  the 
clay  shale.  The  exact  order  for  the  different  materials  changes 
somewhat  depending  upon  which  crater  dimension  is  being  con¬ 
sidered.  A  second  important  fact  that  can  be  noted  is  that  in 
almost  every  case,  nuclear  events  produced  smaller  scaled 
craters  than  their  high-explosive  counterparts  in  the  sam^ 
material.  If  attention  is  brought  to  bear  on  the  surface  burst 
data  for  the  Air  Vent  Series,  which  also  included  the  two  20 
ton  Flat  Top  events,  it  becomes  obvious  that  a  larger  scaling 
factor  would  reduce  the  data  scatter  considerably.  When  this 
fact  is  considered  along  with  a  close  look  at  other  surface 
burst  data,  it  would  appear  that  surface  and  near  surface  bursts 
should  be  studied  separately  and  should  possibly  have  been 
eliminated  from  the  data  used  in  this  research. 

Discussion  of  Approach  Used.  S’ nee  a  good  portion  of  the 
time  spent  on  this  research  was  in  regression  analysis,  it  seems 
appropriate  to  discuss  this  process  and  its  applicability  to 
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research  of  this  type.  In  general,  it  is  sufficient  to  say  that 
multiple  regression  analysis  for  such  a  complicated  problem  as 
was  encountered  here  is  more  of  an  art  than  a  scientific 
approach.  For  regression  analysis  to  be  of  real  benefit, 
experience  with  its  application  and  at  least  some  knowledge  of 
how  the  dependent  variable  behaves  as  a  function  of  the  various 
independent  variables  is  hignly  desired.  Otherwise,  it  becomes 
the  problem  instead  of  a  tool  to  help  solve  the  real  problem. 
Statistics  books,  in  general,  tend  to  oversimplify  the  process 
(rightfully  so  if  they  are  to  teach  the  principles  involved)  but 
the  real  world  just  does  not  seem  to  be  composed  of  only  one  or 
two  independent  variables.  If  a  general  series  is  assumed  as 
the  regression  model,  a  very  few  variables,  applied  to  say 
a  fourth  order  fit,  produces  an  excessive  number  of  terms  to  be 
evaluated.  This  number  of  terms  can  quickly  exceed 
the  number  of  the  observations  or  make  it  extremely  difficult  to 
obtain  an  accurate  solution  of  the  normal  equations  even  using 
the  computer.  It  is  only  when  the  behavior  can  be  simply 
explained  or  when  reasonably  precise  laws  govern  the  behavior 
that  regression  analysis  can  be  used  with  the  assurance  of 
gratifying  results. 

Even  though  the  idea  of  using  the  bell  curve  as  the 
regression  model  was  the  key  to  any  success  achieved  in  this 
research,  it  did  not,  of  course,  fit  all  the  data  well.  Several 
methods  were  tried  in  an  attempt  to  skew  this  model  to  improve 
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the  data  fit,  however  each  method  had  its  inherent  problems. 
Although  the  method  chosen  works  satisfactorily,  it  no  longer 
is  a  bell  curve.  It  is  asymptotic  to  the  x-axis  on  only  one 
end.  If  this  end  falls  toward  the  deeper  depth  of  burst  portion 
of  the  graph  then  no  real  problem  exists  provided  the  "S" 
portion  of  the  curve  does  not  situate  itself  in  the  middle  of 
the  data  as  can  be  seen  for  the  clay  shale  plots  in  Appendix 
III.  Extreme  care,  therefore,  must  be  taken  when  using  this 
model  to  insure  that  the  proper  portion  of  the  curve,  did  in 
fact,  fit  the  data  points  provided. 

As  can  be  noted  from  the  results  of  the  surface  fits,  the 
standard  bell  curve  fitted  more  of  the  data  as  well  or  better 
than  the  skewed  form.  It  also,  in  general,  produced  a  better 
fit  for  scaled  radius  data  than  it  did  for  scaled  depth  and 
scaled  cube  root  of  volume  data.  If  the  assumption  is  made 
that  the  bell  curve  is  truly  representative  of  true  crater 
dimensions,  then  this  fact  makes  sense.  The  apparent  crater 
radius  is  essentially  equal  to  the  true  radius  from  depths  of 
burst  ranging  from  zero  to  past  optimum.  It  is  only  at  the 
deeper  depths  of  burst  that  these  two  radii  diverge.  Of  the 
three  crater  dimensions,  poorer  fits  were  obtained  for  apparent 
crater  depth.  Fall  back  is  the  important  factor  here,  the 
volume  (and  in  turn  the  height)  of  the  fall  back  being  dependent 
on  bulking  and  other  material  characteristics.  Judging  from 
the  data  scatter,  material  properties  are  very  critical  for  this 
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particular  crater  dimension.  Judging  also  from  the  data  plots, 
there  exists  a  range  of  deep  depths  of  burst  in  the  explosion 
process  where  the  material  either  responds  or  does  not  respond 
too  well  to  being  thrown  out. 

Obviously,  from  the  above  discussion,  material  properties 
are  important  in  determination  of  crater  size.  Their  inclusion 
into  the  bell  curve  prediction  scheme  is  covered  next. 


SECTION  VII 


MATERIAL  EFFECTS 

This  section  describes  the  method  used  to  incorporate  soil 
and  rock  properties  into  the  bell  curve  regression  models,  the 
resulting  general  equations  obtained  and  the  analysis  of  these 
equations  to  determine  specific  material  property  effects.  It 
also  includes  a  discussion  of  material  properties  in  general 
and  of  the  material  properties  used  in  this  research  in  partic¬ 
ular.  Lastly,  it  presents  a  brief  discussion  of  the  relation¬ 
ship  between  material  properties  and  cratering  mechanics. 

Incorporation  of  Material  Properties.  After  normalizing 
the  crater  dimensions  and  after  determining  that  the  bell  curve 
was  an  appropriate  regression  model  to  use  if  material  property 
effects  were  ignored,  it  was  surmised  that  the  position  of  the 
bell  curve  (i.e.,  the  height,  point  of  zero  slope,  etc.)  was  a 
function  of  the  media  material  properties.  In  other  words ,  the 
bell  curve  regression  coefficients  obtained  for  the  various 
materials  are  really  a  function  of  appropriate  material  prop¬ 
erties.  This  led  to  regression  analysis  involving  inclusion  of 
soil  and  rock  properties. 

Ideally,  if  sufficient  data  had  been  available  at  various 
material  property  conditions,  then  each  coefficient  (B  value)  of 
the  bell  curve  could  have  been  expressed  as  a  function  of  these 
conditions.  However,  to  make  full  use  of  the  data,  it  was 
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necessary  to  assume  that  the  coefficients  of  the  coded  bell 
curve,  the  A  values,  were  functions  of  material  properties.  The 
B  values  are,  of  course,  related  to  the  A  values  as  described 
by  equations  (5)  and  (6),  but  when  the  B  values  are  shown  in 
terms  of  material  parameters,  the  resulting  equation  is  so 
involved  that  it  becomes  meaningless.  It  was  quickly  determined 
that  the  coded  bell  curve  constants  were  not  linear  functions  of 
any  one  or  more  material  properties.  Reasonable  regression  began 
to  occur  when  second  order  and  interaction  terms  were  included. 
This,  however,  reduced  the  number  of  material  properties  which 
could  be  considered  at  one  time  without  increasing  the  number  of 
parameters  being  used.  Every  effort  was  made  throughout  the 
research  to  keep  the  number  of  parameters  to  40  or  less.  In 
the  final  analysis,  the  best  prediction  formulas  were  obtained 
using  either  equation  (5)  or  (6)  and  the  following  functional 
relationships  for  their  coefficients: 


A-,  =  C-,  +  C2  y5/16  +  C3  S  +  C4  M  +  C5  y5/8  +  C6  Ev  (7a) 
+  C?  M2  +  C8  y5/16  S  +  Cg  y5/16  M  +  C]0  SM 


+  C12 

y5/16  ♦  .  .  . 

•  +  c20  SM 

(7b) 

+  c22 

y5/16  *  .  .  . 

.  +  C30  SM 

(7c) 

+  C32 

V5/’6  ♦  .  .  . 

•  +  C40  $M 

(7d) 

in  which  y  =  total  unit  weight  of  the  material,  in  grams  per 
cubic  centimeter;  S  =  degree  of  saturation,  ranging  from  zero 
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to  1.0;  Ey  =  vaporization  energy  of  the  material,  in  thousands 
of  pounds  per  square  inch  per  cubic  inch  (this  factor  is  zero 
for  all  except  the  nuclear  events);  and  where  M  may  equal  any 

one  of  the  following:  (1)  G$ ,  the  grain  bulk  specific  gravity; 

1  /3 

(2)  tan  <j>,  the  material 's  shearing  resistance;  or  (3)  c  ,  the 
seismic  velocity  of  the  material,  in  feet  per  second. 

Appendix  IV  contains  the  equation  coefficients  (C  values) 
for  the  various  crater  dimensions  ar,d  combinations  of  material 
properties.  Numerous  computer  runs  were  made  to  determine  which 
material  properties  correlated  best  with  crater  dimensions  and 
produced  the  best  surface  fit.  As  could  be  expected,  those 
material  properties  which  were  predominately  measured  values 
correlated  che  best  and  these  included  the  dry  unit  weight  and 
moisture  content  and  to  a  lesser  degree  the  shearing  resistance. 
Dimensional  analyses  by  Westine  (88)  and  Saxe  (61)  suggested 
the  use  of  and  c^3.  These  scaled  values  produced  better 

regression  than  when  their  full  values  were  used.  Degree  of 
saturation,  S,  and  specific  gravity,  G$ ,  were  selected  after 
trying  other  related  factors  such  as  porosity,  percent  air  and 
void  ratio.  They  appeared  to  correlate  better  and  have  a  little 
more  meaning  when  considering  their  effect  on  crater  dimensions. 

The  vaporization  energy  of  the  material,  E  ,  was  used 
primarily  to  differentiate  betv/een  nuclear  and  high  exDlosive 
events.  From  actual  experience  (45)  it  was  known  that  nuclear 
events  produced  smaller  craters  than  would  be  predicted  by 
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scaling  from  high  explosive  cratering  events.  This  fact  became 
evident  when  the  nuclear  events  were  analyzed  with  and  without 
the  £y  term  as  a  parameter.  To  keep  the  number  of  parameters 

2 

to  a  minimum,  Ey  was  used  in  place  of  the  normally  expected  S 
term  in  the  second  order  material  property  scheme  utilized.  The 
S  term  was  found  to  have  the  least  effect  of  all  the  nine 
material  property  terms  on  the  scaled  crater  dimension. 

Also  a;'  attempt  was  made  to  incorporate  all  the  material 
parameters  used  into  one  40  parameter  equation,  by  considering 
the  terms  in  previous  surface  fits  which  had  the  least  effect 
on  the  dependent  variable.  This  attempt,  however,  did  not 
produce  as  good  a  surface  fit  as  when  only  four  material  prop¬ 
erties  were  used. 

When  this  research  was  first  undertaken,  it  was  hoped  that 
crater  dimensions  for  at  least  90  percent  of  the  events  could 
be  predicted  within  ±  10  percent.  After  viewing  the  accuracy 
of  the  basic  cataloged  data,  a  goal  of  80  percent  of  the  events 
to  be  predicted  within  t  20  percent  was  established.  This 
second  goal  was  more  than  met  for  radius  and  volume  but  not 
quite  met  for  predicting  depth.  The  final  figures  obtained  do 
give,  however,  some  indication  as  to  the  reliance  of  using 
material  properties  to  predict  crater  dimensions  and  provide 
some  validity  to  the  assumption  that  investigation  of  the 
general  equations  obtained  would  allow  a  determination  of  the 
effect  of  soil  and  rock  properties  on  crater  dimensions. 
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Although  an  attempt  was  made  to  go  back  over  the  data  to  try 
and  find  specific  explanations  as  to  why  specific  events  did 
not  predict  well,  very  little  success  was  obtained.  In  general, 
it  appeared  that  the  accuracy  of  material  property  inputs  was 
the  predominate  reason.  Not  using  a  better  scaling  factor  for 
the  surface  bursts  accounted  for  a  small  amount  of  discrepancy. 
Even  Palanquin  (a  nuclear  test)  (82)  which  did  not  produce  the 
crater  expected  and  which  in  turn  was  blamed  on  a  stemming 
failure,  predicted  well.  The  general  equations  using  y5/^,  S, 

G  and  E  predicted  Palanquin's  radius  and  volume  only  eight 

o  V 

and  five  percents  higher  respectively  than  those  which  occurred 
and  actually  predicted  the  depth  5  percent  lower  than  the 
actual . 

How  well  will  the  general  equations  obtained  for  the 
various  crater  dimensions  predict  future  events?  Judging  from 
the  data  observed,  there  is  an  85  percent  chance  of  being  within 
±  15  percent  for  materials  that  fall  within  the  range  of  the 
ones  used  in  this  research  and  for  a  scaled  depth  of  burst  less 
than  2. 

Effects  of  Material  Properties.  After  development  of  the 
general  equations  for  crater  dimensions,  parametric  studies  of 
these  equations  were  made  to  determine  how  crater  geometry  was 
effected  as  the  various  material  properties  were  varied  over 
their  normally  expected  ranges.  This  was  accomplished  by 
writing  a  second,  but  much  smaller  computer  program,  to  vary 
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the  parameters  in  the  final  prediction  equations  to  determine 
their  effect  on  final  crater  geometry.  Basically,  this  program 
read  in  the  constants  obtained  for  the  prediction  equations, 
varied  the  material  parameters  and  scaled  depth  of  burst  within 
their  expected  ranges  and  computed  the  estimated  values  for 
scaled  radius,  depth,  or  cube  root  of  volume.  Since  these 
studies  produced  pages  and  pages  of  computer  output,  only 
representative  results  are  included  here.  These  results  are 
graphically  presented  in  Appendix  V. 

Material  Property  Definition.  Because  the  results  and 
conclusions  regarding  the  material  properties  used  in  this 
research  are  presented  next,  it  seems  appropriate  to  discuss 
just  what  is  a  "material  property?"  The  term  "material  property" 
appears  to  have  many  meanings.  These  meanings,  as  they  apply 
to  engineering  problems,  can  be  catagorized  into  three  areas: 

(1)  Basic,  (2)  material  test >and  (3)  theoretical. 

Basic  or  primitive  material  properties  are  those  that 
relate  mass  and  volume.  Examples  for  a  single  phase  material 
are:  (1)  Density,  (2)  specific  volume  and  (3)  specific  gravity. 
For  a  three  phase  material  mixture,  such  as  soil,  where  mineral 
particles  comprise  the  solid  phase  and  where  water  and  air 
occupy  the  voids  between  the  mineral  particles,  the  basic 
material  properties  relate  not  only  mass  to  volume  for  each 
phase  but  relate  mass  and/or  volume  of  each  phase  to  the  total 
volume  or  total  mass.  Examples  of  such  basic  properties  are: 


(1)  Dry  unit  weight;  (2)  total  unit  weight;  (3)  saturated  unit 
weight;  (4)  moisture  content;  (5)  void  ratio;  (6)  porosity; 

(7)  percent  saturation;  and  (3)  percent  air  voids. 

The  second  major  category  is  material  test  properties. 

They  are  defined  operationally.  The  results  of  specified  and 
arbitrary  tests  yield  these  properties.  Examples  of  such  prop¬ 
erties  are:  (1)  Unconfined  compressive  strength;  (2)  Atterberg 
limits  for  soils;  (3)  splitting  tensile  strength  for  rock;  (4) 
compressive  strength  of  concrete  at  28  days;  and  (5)  Hveen 
stability  of  bituminous  concrete.  Such  tests  as  these  are 
usually  standardized. 

There  may  be  some  logic  behind  those  test  procedures  that 
in  some  way  model  a  mechanics  problem,  but  in  essence  these 
test  properties  are  empirical  and  arbitrary  and  may  not  relate 
at  all  to  the  mechanical  behavior  of  the  material  in  other 
situations.  If  they  do,  it  is  indeed  fortunate  and  a  tribute 
to  some  investigator's  intuitive  insight  into  material  behavior. 

The  third  type  of  material  properties  depend  on  a  theory 
of  material  behavior  that  relates  cause  to  effect.  Usually  the 
cause  is  stress,  or  possibly  temperature,  and  the  effect,  strain 
or  rate  of  strain.  The  coefficients  in  these  theoretical 
relationships  are  the  constitutive  constants  that  describe 
idealized  material  behavior.  The  mathematical  model  used  to 
describe  a  material  may  be  as  complicated  as  necessary  to  cover 
the  range  of  interest  of  the  material's  behavior.  Examples  of 


very  simple  constitutive  equations  are  those  for  rigid  bodies, 
perfect  fluids,  linear  elastic  solids,  perfect  gases,  and 
linear  viscous  fluids.  The  constants  in  these  equations  are 
called  material  properties.  For  example,  in  linear  elastic 
theory,  there  are  two  material  properties:  the  two  Lame' 
constants.  More  complicated  constitutive  equations  of  nonlinear 
form  that  relate  the  stress  tensor  to  the  strain  tensor  and 
rate  of  strain  tensor  have  been  generated.  In  the  case  of 
plasticity  theory,  the  constants  relate  the  various  invariants 
of  the  stress  tensor  at  failure.  In  each  of  these  cases,  the 
theory  must  exist  before  a  material  property  is  measured. 

Without  the  theory,  the  property  does  not  exist. 

It  would  be  ideal  and  very  fortunate  if  valid  relationships 
existed  between  the  various  categories  of  material  properties, 
however,  at  best,  we  are  fortunate  that  properties  in  one 
category  may  be  indicative  of  properties  in  another. 

Material  Properties  Used  in  this  Research.  The  best 
regression  equations  were  obtained  as  a  result  of  using  the 
total  unit  weight,  the  percent  saturation,  the  specific  gravity 
of  the  grains,  the  internal  ..hearing  resistance  and  the  seismic 
velocity.  Although  these  parameters  may  not  really  be  the 
properties  governing  material  behavior  during  cratering,  they 
are  at  least  indicators  of  the  properties.  None  of  these 
parameters  except  the  specific  gravity  remain  constant  during 
the  cratering  process.  They  are  all  functions  of  the  stress 
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(or  strain)  level  and  their  initial  values  are  not  necessarily 
indicative  of  their  values  during  the  actual  explosive  event. 
However,  when  used,  in  conjunction  with  one  another  and  in 
conjunction  with  the  depth  of  burst,  they  provide  some  indication 
of  the  governing  material  properties.  It  was  interesting,  how¬ 
ever,  to  study  tie  effect  these  indicators  did  have  on  crater 
dimensions . 

At  the  optimum  depth  of  burst  for  granular  materials , 
larger  craters  were  produced  at  low  and  high  degress  of  satu¬ 
ration  than  were  produced  at  the  intermediate  values  (Figs.  32- 
34).  This  was  probably  because  granular  materials 
with  low  moisture  contents  have  very  little  cohesion.  However, 
as  the  moisture  content  is  increased,  the  material  develops 
some  cohesion  and  therefore  greater  strength.  Whitman  (90) 
showed  that  cohesion  due  to  pore  water  increased  as  the  rate 
of  strain  increased.  The  addition  of  water  also  increased  the 
weight  of  the  material.  When  the  degree  of  saturation  starts 
to  approach  100  percent,  the  ability  of  the  material  to  transmit 
the  shock  wave  markedly  increases.  In  addition,  the  strength 
of  the  material  decreases  considerably  as  the  pore  water  starts 
to  assume  more  and  more  of  the  pressure  being  exerted  on  the 
mate.ial.  For  clays  and  rocks  (Fig.  35),  any  increase  in  the 
degree  of  saturation  appeared  to  produce  larger  craters. 
Additional  moisture  in  these  cases  did  not  improve  their 
cohesion  but  oniy  tended  to  increase  their  ability  to  transmit 


the  shock  wave  and  to  decrease  their  strength. 

Now,  if  other  than  optimum  depth  of  burst  is  considered, 
the  degree  of  saturation  produced  slightly  different  effects 
(Fig.  32).  For  deep  depths  of  burst,  an  increase  in  the  degree 
of  saturation  tends  to  always  produce  larger  craters.  For 
surface  bursts,  the  effect  of  the  degree  of  saturation  was 
found  to  be  the  reverse  of  that  determined  for  optimum  depth 
of  burst.  Low-  and  high-moisture  contents  produced  smaller 
craters,  and  there  existed  an  optimum  degree  of  saturation  at 
which  the  largest  crater  was  obtained.  Since  compaction  is  the 
predominate  mechanism  here,  it  follows  that  materials  in  this 
region  would  behave  essentially  as  they  do  in  normal  engineering 
compaction  problems. 

Again,  considering  cratering  events  at  or  near  optimum 
depth  of  burst  (Figs.  36  and  37),  it  was  found  that  the  smaller 
craters  resulted  from  low  density  and  high  density  materials 
and  that  there  existed  an  optimum  density  at  which  the  largest 
crater  resulted.  This  feature  reflects  again  that  better  prop¬ 
erties  are  necessary  to  predict  crater  dimensions.  If  unit 
mass,  strength  and  the  ability  of  the  material  to  absorb  the 
shock  wave  are  all  considered  together,  then  this  phenomenon 
seems  plausible.  Materials  with  low  weight  and  strength  appear 
to  have  high-energy  absorption  properties  whereas  dense,  high- 
strength  materials  do  not.  Somewhere  in  the  middle,  then, 
there  exists  a  material  where  these  factors  do  not  compliment 
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each  other  as  much  as  they  do  for  the  two  extremes.  This  density 
feature  appears  to  hold  true  for  all  depths  of  bursts. 

Because  the  effect  of  bulk  grain  specific  gravity  on  crater 
dimensions  is  not  so  meaningful  in  terms  of  its  application  to 
the  effects  of  material  properties,  it  will  only  be  discussed 
briefly.  A  change  in  grain  specific  gravity  affects  the  total 
unit  weight  and  degree  of  saturation  and  those  were  discussed 
above. 

Not  so  easily  explained  as  the  effect  of  density  and  degree 
of  saturation  on  cratering  is  the  effect  of  the  material's 
internal  shearing  resistance,  tan  <j>.  It  is  very  difficult  to 
determine  what  the  actual  ranges  of  tan  <f>  are  for  a  particular 
soil  density  and  degree  of  saturation.  Every  effort,  then,  was 
made  to  stay  within  the  area  of  the  actual  data  to  determine 
the  effect  of  tan  <p  on  crater  dimensions.  In  general,  it  appeared 
that  as  tan  <p  increased  for  soil  materials  at  the  lower  degrees 
of  saturation  and  for  the  rocks,  the  size  of  the  resultant 
crater  also  increased.  If  high  degrees  of  saturation  are  con¬ 
sidered  for  the  soils,  then  an  increase  in  tan  <j>  decreases  the 
crater  size.  The  only  explanation  that  seems  plausible  for 
these  effects  is  to  consider  the  ability  of  the  material  to 
absorb  energy  in  relation  to  tan  <f>.  Apparently  the  ability  of 
the  material  to  absorb  energy  increases  more  with  an  increase 
in  internal  shearing  resistance  than  does  the  strength  for  the 
lower  degrees  of  saturation.  At  the  high  degrees  of  saturation, 


however,  the  energy  absorption  value  apparently  levels  off  and 
the  strength  increase  is  sufficient  to  produce  smaller  craters. 

Because  the  parametric  study  involving  tan  <j>  appeared  to  be 
somewhat  meaningless  unless  values  were  considered  in  the  same 
area  as  the  actual  data,  a  parametric  study  was  not  performed 
using  the  general  equation  which  included  seismic  velocity. 

Again  it  becomes  extremely  difficult  to  determine  the  range  of 
values  in  seismic  velocity  which  would  exist  when  the  density 
and  degree  of  saturation  are  assumed.  In  addition,  there  were 
fewer  measured  values  for  this  parameter  and  it  was  included  in 
the  general  equations  with  hesitancy.  It  was  felt,  however, 
that  even  though  the  data  were  not  the  best  and  this  variable 
is  stress  dependent,  it  very  likely  would  give  seme  indication 
of  the  material's  ability  to  transmit  the  shock  wave.  A 
cursory  review  of  the  data  indicates  that  smaller  craters 
result  when  the  seismic  velocity  is  either  low  or  high  and  that 
there  exists  an  optimum  value  where  the  largest  crater  will  be 
produced. 

It  becomes  very  difficult  to  sum  up  all  the  possibilities, 
but  in  general  it  appears  that  for  craters  produced  at  or  near 
optimum  depth  of  burst,  there  exists  an  optimum  unit  weight,  an 
optimum  internal  shearing  resistance,  an  optimum  seismic 
velocity  with  the  degree  of  saturation  at  zero  percent  where 
the  very  largest  craters  will  be  produced  for  a  particular 
explosive  energy  source. 


53 


For  practical  considerations,  it  appears  that  the  most 
important  of  the  material  indicators  used  is  the  degree  of 
saturation  (which  is  of  course  the  result  of  water  content). 

This  study  should  be  helpful  in  determining  whether  the  addition 
of  water  to  the  medium  will  either  enhance  or  decrease  the  size 
of  a  crater  being  considered  for  the  explosive  charge  weight 
being  used. 

Cratering  Mechanics  and  Material  Properties.  As  a  result 
of  this  analysis  of  the  effects  of  material  properties  on  crater 
dimensions,  it  appears  very  likely  that  a  reasonably  simple 
cratering  theory  should  be  possible.  Hopefully  the  coefficients 
(material  properties)  in  this  theoretical  constitutive  relation¬ 
ship  could  be  measured  using  simple  laboratory  or  field 
techniques.  In  any  case  it  would  appear  that  these  coffi dents 
should  in  some  way  be  related  to  the  following  material  prop¬ 
erties:  (1)  Energy  dissipation;  (2)  total  unit  weight;  (3)  shear 
strength;  (4)  volume  change;  and  (5)  moisture. 

These  five  material  properties  should  be  measured  over  the 
material  field  for  each  future  test  as  a  minimum  material 
property  requirement.  The  energy  dissipation  constant  should 
account  for  the  fact  that  smaller  craters  result  in  the  light¬ 
weight,  weaker  materials.  Unit  weight  in  conjunction  with  depth 
of  burst  would  give  a  measure  of  the  total  mass  of  the  material 
which  must  overcome  gravity.  The  shear  strength  constant  would 
account  for  further  energy  dissipation.  A  volume  change 
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constant  would  account  for  the  change  in  density  of  the  crater 
fall  back  which,  in  turn,  effects  the  apparent  crater  depth  and 
volume,  Last  but  not  least  is  the  moisture.  Vaporiation  of 
moisture  surrounding  the  explosive  charge  enhances  the  gas 
acceleration  pi^se  of  cratering  mechanics.  In  addition, 
moisture  would  modify  the  effects  of  all  the  other  four 
properties  proposed. 


SECTION  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 

If  explosives  are  to  be  used  for  excavation  purposes,  then 
predictions  of  results  must  incorporate  material  properties. 

This  study  has  shown  that  soil  and  rock  properties  are  important 
ir.  determining  the  size  of  explosion-produced  craters  and  has 
provided  some  insight  as  to  their  specific  effects.  It  has 
shown  that  previous  investigators  have  been  somewhat  negligent 
in  measuring  material  properties  fcr  past  cratering  experiments 
and  that  no  real  analysis  of  crater  data  can  be  made  unless  the 
variables  are  either  controlled  or  measured.  This  study  also 
has  provided  a  means  to  predict  crater  dimensions  for  any 
material  provided  certain  soil  and  rock  properties  are  measured 
beforehand. 

The  final  general  equations  obtained  will  predict  the  size 
of  85  percent  of  the  cratering  events  within  ±  15  percent. 

These  equations  predict  apparent  crater  radius,  depth,  and  volume 
in  terms  of  (1)  depth  of  burst,  (2)  the  explosive  weight  and 
(3)  the  following  material  properties:  total  unit  weight, 
degree  of  saturation,  grain  specific  gravity,  internal  shearing 
resistance  and  seismic  velocity.  For  nuclear  explosions,  the 
vaporization  energy  of  the  material  is  also  included.  These 
equations  were  developed  from  data  which  fell  in  the  following 
ranges:  (1)  Explosive  charge  weights  from  one  pound  to  100 
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kilotons;  (2)  materials  in  which  the  unit  weight  ranged  from 
60  to  170  pounds/cubic  foot;  and  (3)  charge  depths  from  zero 
to  the  point  where  no  crater  is  produced. 

It  is  not  surprising  that  previous  investigators  have 
concluded  that  the  effect  of  material  properties  on  crater 
dimensions  was  somewhat  obscure.  This  is  particularly  true  if 
only  linear  and  simple  curvilinear  relationships  are  considered. 
It  is  only  when  indicative  soil  and  rock  properties  are  con¬ 
sidered  in  conjunction  with  one  another  for  specific  ranges  of 
depths  of  burst  that  they  become  meaningful. 

Of  the  six  material  properties  used  in  the  general 
equations,  percent  saturation  and  total  unit  weight  appeared 
to  be  the  most  important  indicators  of  the  effects  of  material 
properties.  This  was  due  primarily  because  these  two  properties 
were  calculated  from  predominately  measured  values  while  those 
for  the  other  properties  used  were  largely  estimated.  Prime 
examples  of  these  effects  are  as  follows:  (1)  For  surface 
explosions,  there  exists  an  optimum  percent  saturation  which 
will  produce  the  largest  crater;  (2)  for  soils  and  for  the 
explosive  charge  at  optimum  depth,  zero  moisture  content  pro¬ 
duces  the  largest  crater;  in  this  case  there  is  a  least  favor¬ 
able  percent  saturation  which  will  produce  the  smallest  crater; 
(3)  for  rock  materials  and  for  all  materials  at  deep  depths  of 
burst  (at  least  twice  the  optimum),  100  percent  saturation 
results  in  the  largest  crater,and  (4)  in  general,  there  exists 
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an  optimum  unit  weight  which  will  produce  the  largest  crater. 

The  bell  curve  provides  a  good  model  for  predicting  scaled 
crater  dimensions  in  terms  of  scaled  depth  of  burst.  Where 
data  indicate  nonsymmetry,  the  skewed  form  of  the  bell  curve 
can  be  used,  provided  it  is  done  cautiously.  The  inherent 
advantage  to  the  bell  curve  model  is  that  the  constants  in  this 
model  immediately  provide  the  maximum  value  for  the  scaled 
crater  dimension  under  consideration  and  the  optimum  scaled 
depth  of  burst  at  which  it  occurs.  This  feature  makes  it  useful 
for  practical  applications  and  allows  quick  comparisons  between 
rock  and  soil  property  conditions. 

As  a  result  of  the  knowledge  gained  in  this  research 
effort,  the  following  recommendations  are  made  for  future 
cratering  experiments  and  associated  studies: 

1.  Sufficient  material  properties  should  be  measured 
for  each  and  every  cratering  event.  As  a  minimum, 
unit  weight,  moisture  content,  grain  specific 
gravity,  shearing  resistance  and  seismic  velocity 
should  be  measured.  In  addition,  some  measure  of 
energy  dissipation  and  material  bulking  would  be 
desirable. 

2.  A  method  needs  to  be  developed  to  measure  the 
energy  dissipation  characteristics  for  soil  and 
rock  in  which  cratering  experiments  have  and  will 
be  performed. 
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3.  A  method  to  measure  and  evaluate  bulking  or 
compaction  needs  to  be  developed. 

4.  A  series  of  laboratory  type  experiments  to 
consider  specifically  the  effects  of  percent 
saturation  on  crater  dimensions  would  be  highly 
desirable.  As  a  minimum,  five  levels  of  the 
degree  of  saturation  for  five  levels  of  depth  of 
burst  for  a  particular  dry  unit  weight  would 
possibly  suggest  a  more  accruate  relationship 
between  crater  dimensions  and  this  very  important 
material  parameter. 

This  series  could  also  be  extended  to  include 
a  variation  in  the  dry  unit  weight  and  the 
inclusion  of  various  types  of  materials. 

5.  For  survival  of  silo-launched  missile  systems, 
where  the  primary  interest  is  in  near  surface 
bursts,  the  surface  burst  data  from  this  study 
could  be  supplemented  with  additional  date  and  a 
regression  analysis  performed.  This  would  allow 
the  scaling  exponent,  as  well  as  the  material 
constant,  to  be  a  function  of  material  properties. 
Thus,  a  more  accurate  prediction  of  crater 
dimensions  could  be  obtained  for  this  special  case. 

6.  A  simple  theory  of  cratering  should  be  developed 
using  the  field  equations  of  mechanics  and  a  material 
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constitutive  equation  with  sufficient  complexity 
to  account  for  soil  or  rock  failure  and  energy 
absorption.  Such  a  constitutive  equation  should 
also  account  for  increase  in  strength  as  a  functi 
of  pressure  and  the  thermodynamic  properties. 


APPENDIX  I 


CATALOGED  CRATER  DATA 

Notations  and  Definitions.  The  following  notations  and 
abbreviations  are  used  in  the  cataloged  crater  data  and  material 
property  data  listings: 

ALLUV  -  desert  alluvium; 

ATTBRG  LIMITS  -  Atterberg  limits  (see  also  LL  and 
PI)  relating  to  the  water  content  at  which 
soil  consistency  changes  from  one  state  to 
another,  see  Wu  (91); 

B  -  indicates  that  there  is  no  value  following  even 
though  the  computer  printed  zeros; 

BULK  FACTOR  -  bulking  factor,  a  ratio  of  the  unit 
weight  of  the  material  in  the  crater  fallback 
to  the  preshot  unit  weight  of  the  material. 

Frandson  (19)  analyzed  this  value  for  several 
materials; 

CH  -  inorganic  clays  of  high  plasticity; 

CNF  PRES,  PS  I  -  confining  pressure  at  which  the 
confined  compressive  strength  was  obtained,  in 
pounds  per  square  inch; 

COHESION,  PSI  -  cohesion  of  the  material,  in  pounds 
per  square  inch,  based  on  Mohr-Coulomb  failure 
theory; 


CONF  COMP,  PSI  -  confined  compressive  strength  of 
the  material,  in  pounds  per  square  inch,  at 
the  particular  confining  pressure  listed  ir 
the  next  column; 

CORE  RECOV,  PRCT  -  the  amount  of  core  recovered 
during  coring  operations,  in  percent;  for  rock 
it  indicates  the  material's  soundness; 

CU  FT  -  cubic  feet; 

DRY  UWT,  L8S/CU  FT  -  dry  unit  weight  of  the  material 
in  pounds  per  cubic  foot; 

ELEV,  FT  -  elevation,  in  feet,  which  could  have  been 
converted  to  atmospheric  pressures  which  Herr 
(25)  showed  to  be  important; 

EQUIV  WT,  LBS-TNT  -  equivalent  of  the  explosive 
charge  in  pounds  of  trinitrotoluene; 

EVT  NO.  -  event  number; 

FL.E  STN,  IN/IN  -  the  strain  at  which  failure  occurred 
in  either  the  unconfined  or  confined  compression 
test,  in  inches  per  inch; 

FPR  MNT  -  Fort  Peck  Reservoir,  Montana; 

FT  -  feet; 

KT  -  kiloton,  one  thousand  tons  equivalent  weight 
of  TNT  (trinitrotoluene); 

LIP  HT,  FT  -  apparent  crater  lip  height,  in  feet; 
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LL,  PR  CT  -  liquid  limit  (see  ATTBRG  LIMITS) ,  in 
percent;  the  water  content  of  the  soil  which 
differentiates  between  the  plastic  and  liquid 
consistencies  of  the  soil; 

LRL  -  Lawrence  Radiation  Laboratory,  California; 

MELT,  MPSI/CIN  -  the  energy  required  to  melt  the 
material,  in  thousands  of  pounds  per  square  inch 
per  cubic  inch; 

ML  -  inorganic  silts  and  very  fine  sands,  silty  or 
clayey  fine  sands  or  clayey  silts  with  slight 
plasticity; 

MOIST,  PR  CT  -  moisture  content,  expressed  in  per¬ 
cent  of  the  dry  unit  weight  of  the  material; 

MTCE  -  Multiple  Threat  Cratering  Experiment; 

NM  -  nitromethane; 

NTS-A5  -  Nevada  Test  Site  Area  5; 

NUC  -  nuclear; 

PHI,  DEG  -  the  angle  of  internal  shearing  resistance, 
in  degrees  , based  on  Mohr-Coulomb  failure  theory; 

PI,  PR  CT  -  plasticity  index  (see  ATTBRG  LIMITS),  in 
percent,  the  water  content  difference  between  the 
liquid  limit  (LL)  and  the  plastic  limit;  the 
plastic  limit  being  that  water  content  of  the  soil 
which  differentiates  between  the  semi  sol  id  and 
and  plastic  consistencies  of  the  soil; 
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POISN  RATIO  -  Poisson's  ratio,  ratio  of  the  horizontal 
stress  to  the  vertical  stress,  which  resulted  from 
the  theory  of  elasticity; 

PRE-GOLA  -  Pre-Gondola; 

REF  NO.  -  reference  number; 

RMK,  SEE  NTE  -  remarks,  see  note; 

SLP  DEG  -  approximate  angle  the  apparent  crater  slope 
makes  with  the  horizontal  preshot  ground  surface; 

SM  -  silty  sands,  sand-silt  mixtures; 

SP  -  poorly  graded  sands,  gravelly  sands,  little  or 
no  fines; 

SP  GR  -  bulk  specific  gravity  of  soil  or  rock  grains; 

TNS  -  tons; 

TNSLE-D,  PSI  -  direct  tensile  strength  of  the  material, 
in  pounds  per  square  inch; 

TNSLE-S,  PSI  -  splitting  tensile  strength  of  the 
material,  in  pounds  per  square  inch; 

UNC  COMP,  PSI  -  unconfined  compressive  strength  of 
the  material,  in  pounds  per  square  inch; 

USCS  CLASS  -  the  "Unified  Soil  Classification  System" 
(73); 

VAPOR,  MPSI/CIN  -  the  energy  required  to  vaporize  the 
material,  in  thousands  of  pounds  per  square  inch 
per  cubic  inch; 
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WT -VOL  -  weight-volume; 

YFC  WSH  -  Yakima  Firing  Center,  Washington; 

$  -  indicates  the  value  following  is  an  estimated 
value. 

Cataloged  Crater  Data.  The  crater  data  cataloged  is 
computer  listed  in  Table  2. 

Cataloged  Material  Property  Data.  The  two  line  computer 
listing  of  the  cataloged  material  property  data  associated  with 
the  cataloged  crater  data  list  is  presented  in  Table  3. 

Notes .  The  notes,  referred  to  in  either  the  catalog  of 
crater  data  or  the  catalog  of  associated  material  properties, 
follow  Table  3. 
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TABLE  2.  CATALOGED  CRATER  DATA  (CONTINUED) 

IDENTIFICATION  EXPLOSIVE  OATA  OEPTH  APPARENT  CRATFR  DIMENSIONS  RNX 


76  23  PRE-GOLA  I-CHAR  FPR  HNT  OCT66  SHALE  NM  19.62  TNS  <.3160.00  42.49  80.40  >2.60  277550.00  29  14.50 


TABLE  3.  CATALOGED  MATERIAL  PROPERTY  DATA  (CONTINUED) 


0002 


TABLE  3.  CATALOGED  MATERIAL  PROPERTY  DATA  (CONTINUED) 


0003  J  OOOC  S  000/5  S  0000 


TABLE  3.  CATALOGED  MATERIAL  PROPERTY  DATA  (CONTINUED) 


?oo  o  o.o  b  o.o  r  ioo 


00  T  I  J  0081  1  OCOOil  1  OOOOlZT 


TABLE  3.  CATALOGED  MATERIAL  PROPERTY  DATA  (CONTINUED) 

MATERIAL  MT-VOL  RELATIONSHIPS  STRENGTHS  AN1  STRAINS 

UM  USCS  R£F  SP  GR  ORT  UWT  HOIST  TNSLE-S  TnSLE-0  UNC  COHP  FLE  ST*i  CONF  COMP  CNF  PRES  FLE  STN 


2300 


00  f  2 


APPENDIX  II 


THE  COMPUTER  PROGRAM 

Description  of  the  Program  Elements.  Considerable  time 
and  effort  was  expended  in  developing  the  program  and  its 
formats.  Therefore  it  seems  appropriate  to  include  it  as  an 
appendix.  A  brief  description  of  its  essential  features 
follows. 

Main  Program.  -This  portion  of  the  program  is  nothing  more 
than  a  calling  program,  i.e.,  it  calls  the  various  main  sub¬ 
routines  to  read,  sort  and  print  the  data  and  to  perform  the 
surface  fit  of  selected  data. 

Subroutine  REDATA.  -This  subroutine  reads  and  stores  all 
crater  and  material  property  data  into  the  computer  data  banks. 
Six  cards  are  read  for  each  cratering  event.  In  addition,  any 
note  cards  associated  with  an  event  are  also  read  and  stored. 

Subroutine  PRDATA  and  Related  HEAD  Subroutines.  -This  sub¬ 
routine  prints  all  crater  and  material  property  data  along  with 
calling  the  related  CDHEAD,  MPHEAD1  and  MPHEAD2  subroutines  to 
provide  the  necessary  headings  for  the  output  data.  This  sub¬ 
routine  produces  the  cataloged  data  listing. 

Subroutine  SURFIT.  -This  subroutine  is  the  main  program 
for  the  conduct  of  the  least  squares  surface  fit  and  analysis. 
It  primarily  calls  the  various  subroutines  necessary  to  perform 
the  surface  fit  and  analysis.  It  also  specifies  the  dependent 
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variable  (radius,  depth  or  volume),  the  scaling  exponent  and 
the  number  of  independent  variables  to  be  used  in  the  surface 
fitting  process.  In  addition,  it  calls  for  or  specifies 
printing  of  all  pertinent  matrices  and  coefficients. 

Subroutine  EQN.  -This  subroutine  determines  the  type 
regression  and  the  independent  variables  to  be  used  in  the  least 
squares  surface  fit.  It  develops  the  vector  of  observations  and 
the  matrix  of  measured  independent  variables  and  normalizes 
these  matrices  to  provide  better  matrix  inversion  and  manipu¬ 
lation. 

Subroutine  COP.CO.  -  This  subroutine  develops  a  simple 
correlation  matrix  between  all  the  variables  involved  in  the 
surface  fit;  dependent  as  well  as  independent. 

Subroutine  COEFF .  -This  subroutine  takes  the  vector  of 
observations  and  the  matrix  of  measured  independent  variables 
and  generates  the  vector  and  x-coefficient  matrices.  These 
latter  matrices  are, in  essence,  the  normal  equations  which  must 
be  solved  simultaneously  to  obtain  the  regression  coefficients. 

Subroutine  ABPRN.  -This  subroutine  prints  the  coefficients 
of  the  prediction  equation.  It  also  prints  the  decoded  form  of 
these  coefficients  when  the  model  used  was  either  a  three  or 
four  parameter  bell  curve. 

Subroutine  PREDIC.-This  subroutine  calculates  the 
estimated  value  of  the  dependent  variable  for  each  observation 
using  the  empirical  equation  generated  and  compares  this  value 
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with  the  actual.  It  also  calculates  the  multiple  correlation 
coefficient  and  standard  deviation  for  all  the  data  being 
considered. 

Subroutine  MATINV.  -This  subroutine  inverts  the  /-coeffi¬ 
cient  matrix  to  be  used  in  solving  the  normal  equations. 

Subroutine  MATPRN.  -  This  subroutine  prints  a  square  matrix. 
It  is  used  to  print  the  x-coefficient  matrix  as  well  as  the  unit 
matrix  which  should  result  when  the  x-coefficient  matrix  and  the 
inverse  of  the  x-coefficient  matrix  are  multiplied  together. 

Subroutines  SIMALT  and  MULT.  -These  subroutines  multiply  a 
square  matrix  times  a  column  matrix  and  a  square  matrix  times 
another  square  matrix  respectively.  They  are  used  to  multiply 
the  inverse  of  the  x-coefficient  matrix  times  the  vector  matrix 
to  obtain  the  coefficients  of  the  curve  fit  and  to  multiply  the 
x-coefficient  matrix  times  its  inverse. 

Typical  Program.-  A  computer  printout  of  the  program  as  it 
wa s  run  during  the  latter  stages  of  the  research  follows. 

Sample  Data  Output.-  The  computer  output  for  a  simple 
regression  analysis  of  scaled  radius  as  a  function  of  scaled 
depth  of  burst  and  one  material  property,  total  unit  weight, 
using  the  skewed  bell  curve  model  follows  the  program. 
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tsmmsimsssiifsmisimmmmmmmimsmsmmmmissmmm 


PROGRAM  MAI N< INPUT .OUTPUT, T APES = INPUT, TAPE6=0UTPUT , PUNCH, TAPE7= 
♦PUNCH) 

MAIN  PROGRAM 

CALL  REOAIA  (NOA 7 A) 

CALL  SURF I r ( NO A 7 A ) 

STOP 

ENO 


SUOROUTIMi  REOAT  A (NOATA ) 

C  THIS  SU3R0UT INE  REAOS  ALL  CRATER  ANO  MATERIAL  PROPERTY  DATA  FOR  EACH 

C  CRATERING  EVENT.  SIR  CARDS  ARE  READ  FOR  EACH  EVENT ( I T EH  NO.)  ALONG 

C  WITH  ANY  RENARKS(NOTE)  CAROS.  A  OLANK  CARD  MUST  GE  INSERTED  GETWEEN 

C  OATA  CAROS  ANO  NOTE  CAROS  ANO  AT  THE  END  OF  THE  NOTE  CAROS. 

COMMON  WI200  ,28)  .SHOT  (230,4>,SITE(20  0,?.),DATE(200,2), 

»EMEO(200,2|  ,  EX  TYPE  (20  0)  ,  YIELO(  2  00,3  )  ,  TNT  NT  (200!  ,006(200)  , 

♦RADIUS (200) , DEPTH (200) , VOL (200) ,HTL  IP (200)  , 

♦RMKCO (200) , £LEV(200) , CL  ASS ( 200 , 2) 

COMMON  SPGR (  200) ,UWT (200) , PMOIST(2DO) ,SPTEN(203)  ,0ITEN(200)  , 
♦UCGMP(200) ,UCSTN(200) ,CCOHP(200) ,CONFP(200) ,CCSTN(200) , PHI (200) , 
♦COHES(200I  ,  POI  s.N  (200  '  .BULK  (200)  ,  SEC  HOD  (2  00)  ,YOUMOD(200!  . 
•SHEMOD(200)  .SEIVELUOOI  ,SHEVOL(2QO) , AT  ILL (200 1  , 

♦ATTPI (2001 ,CQRE(200) ,£MELT (200) .VAPOR (2 00) , RMKHP ( 200, 2) 

COMMON  I NO ( 200 ) , NREF (200),ISL(200),KREF(200) 

COMMON  /NOTE/  NCAR0I103) « RNTE ( 100) , TEXT (100,19) 

1  =  0 

100  i  =  m 

READ  (5, 10)  INO(I)  , NREF  (I)  ,  (SHOT  (I,  J)  ,  J=1 ,  4)  ,  SITE  ( I  ,  J)  ,  J=1 , 2 )  . 

♦  ( OATE  ( I ,  J)  ,  J  =  l,2)  ,  (EMEO  (I,  J)  ,  J  =  l,2)  ,  EX  TYPE  ( I )  ,  (  YIELD  ( I  ,  J)  ,  J=l ,  J )  , 
♦TNTHT(I) 

IF(IMOCI) .EQ.OIGO  TO  230 

REAO  (5,201006(1)  ,  RAOIU3  ( I )  ,  OEPTKt  I)  ,  VOL  ( I )  ,  W  ( 1 ,  1)  ,  ISL  (  I)  ,  H  ( 1 , 2)  , 
♦HTLIP(I) ,RMKCO(I) 

REAO  (5, 301  4(  1,3)  ,  ELEV  ( I )  ,  ( CLASS  ( I ,  J )  ,  J  =  1, 2  )  ,HRE,(tf  ,4(1,4), 

♦SPGR (I), HI  1,5) , 0  HT ( I )  ,4(1,6) , PHOIST ( I ) 

Kt  AO  ( 5 ,40)  H  (  1 , 7)  ,SPTEN(  I)  ,H(I  ,8)  ,DITEN(  I)  ,H(I,9)  ,  UCOMP  ( I )  ,  H(  1 ,  1 0)  , 
♦UC3TNII) ,H(I ,11) ,CCOMP( I) ,W(I, 12) ,CONFP (I) ,  W<  1 , 13)  ,  CCSTN  ( I ) 
KEAO(5,50)W(I, 14) , PHI (I ) ,W(1, 15) , COHES ( I) , W( I , 16)  ,  POISNt I ) , 
*W(I,17),flJ_K(I)  ,  W  ( I  ,  1  8>  ,SECMOO(  I)  ,4(1, 19)  ,  YOUMOO(I)  ,H(  1,20)  , 
♦SHEMOD ( I ) 

REA  0(5, 60) H( 1,21) ,SEI VEL (I) ,Ht 1,22) , SHE  VOL  (I) ,H(I, 23) , 

»  H  ( 1 ,24)  ,  AT  TLL(I)  ,H(  I,  25)  ,  ATTPI  ( I)  ,H(I,2  6)  ,  CORE  ( I )  ,  H  ( I ,  27 )  , 
•EHELT(I) ,4< 1,28) ,VAPOR( I) , (RMKHP (I. J) , J=l,2) 

10  FORMAT (IX, 2I4,6A4,2A3 , 2 A 4 , 2X , 3 A4, A 2 , F 1 3 . 2 ) 

23  FORMAT (SX, 3F 3 . 2 , FI  3 , 2 , A  2 , 1  3, A2 , F6.2 , A4 > 

3  0  FORHAT(21X, A 2 , =6 . 0 , 8X , 2 A 3 , 14 , A  2 , F5 . 2 , A2 ,F 6  ,  i , A2 , F5 . 1 ) 

40  FORHAT(5X,2(A2,r6.0) , A2,F7.0,A2,F6.4,A2,F8.0,A2,F7.0,A2,F6.4) 

50  FORHAT(5X,A2,F5.1,A2,F2.0,A2,F4.2,A2,F5.2,3(A2,F9.0) ) 

60  FORHAT(5X,2 (A2.F7.0) ,A2,7X,2(A2,F5.i),A2,F4.0.A2.F6.0,A2,F7.0,A4,A 
*3) 

GO  TO  100 

200  NOA  TA  =  I- 1 
J  =  0 

300  J=JU 

REAO (5,70) NCARO(J) , RNTE ( J) , ( TEX T ( J, K) , K= 1, 19) 

70  FORMAT(I1.F4.0,18A4,A3) 

IF(NCAR0(J) .NE.OJGO  TO  300 

RETURN 

EHO 


SUBROUTINE  PROftTA (NOATA) 

C  THIS  SUBROUTINE  PRINTS  ALL  OATA  REAO  INTO  THE  COMPUTER  3Y  THE  RECATA 
C  SUBROUTINE. 

COMHON  W( 200,28) .SHOT (200, A)  ,  S  IT  E  (200 , 2)  ,  DATE  (200, 2)  , 

•EMEO (200,2)  .EXT Y PE (20  0) , YIELD ( 200, J) ,TNTWrt200),DOB(200), 
•RAOIUS1200) , OEPTH(2QO). VOL (200) ,HTLIP(200)  , 

•RHKCO (200 )  , ELEV ( 200 ) , Cl  ASS (200 , 2) 

COMMON  S°GR( 200 ) ,UHT (20  0) , PMOI ST ( 20  0) ,SPTEN(200)  ,OITEN(200) , 
*UCONP(200) , UCSTN (200 ) ,CCOMP(200) ,C3NFP (200) ,CCS t N( 200) ,PHI (20 0) , 
•COMES (200)  ,POISN(200) .BULK (20  0  ) , SEC  HO 0(200)  .YOUKOO (200) , 

•SHEMOO (200 ) , SCI VEL (20 0) ,SH£VOL(200I  ,ATTLL(200), 

•ATT  PI (200) .CORE (200) , SMELT (200) , VAPOR ( 200) , RHKHP (200,2) 

COHMON  INO(200),NREF(200) , ISL ( 200) , HREF (200) 

COMMON  /NOTE/  NC ARD ( l 03 ) ,RNTE ( 100 ) , TEXT ( l 0 0 , LB) 

1  =  0 

LOO  CALL  COHEAO 
NLINES=7 
200  1*1*1 

WRITE  (6, 10 1 INO(I) ,NREF( I) , (SHOT (I, J) ,J=1,4) , (SITEd , J) ,J=1,2)  , 
*(OATE(I,  J)  ,  J=l,2>  ,  (EMEO  (I ,  J)  ,  J  =  1 , 2)  ,  EXT  YpE  ( I)  ,  (  Y  IELD<  I .  J)  ,  J=1 , 3)  , 
•TNTWT (I) ,003 (I) ,RAOIUS( I) ,OEPTH(I) , VOL (I) , W(I,1) , ISL(I) ,H( I ,2) , 
•HTLIP(I) .RMKCO(I) 

IF ( I . EQ. NOA TA) GO  TO  300 
NLINES=NLINES*2 
IF(NLINES.EQ.59)G0  VO  100 
GO  TO  200 
300  1*0 

AOO  CALL  MPHEOl 
NLINES=8 
K*0 

500  1=1*1 
K=K*1 

WRITE  (6,20)  I  NO  (I)  ,  (SHOT  ( I ,  J)  ,  J  =  l,  4)  ,  W  ( 1 , 3  )  ,  ELEV  { I )  , 

• ( EMEO ( I, J) , J  =  1 , 2) , (CLASS(I.J) ,J  =  1.2) ,HREF(  I) , W ( I , A )  ,SPGR(I) , 
*W(I,5),UWT(I),W(I,6> ,PMOIST(I) ,W(I,7) ,SPTEN(t) ,W(1 ,8) , OITEN(I)  , 
»W(I,9) ,UCOHP(I) ,W(I,10l ,UCSTN(I),W(I,11) ,CCOMP(I) , H(I, 12) , 
»CONFP(I),W(I,13) ,CCSTN(I) 

IF ( I.EQ.NOATA)GO  TO  550 
NLINES=NLINES*2 
IF(NLINES.LT.30)GO  TO  500 
550  CALL  HPHE02 

NLINES=NLINES*7 
I*I-K 
560  1=1*1 

WRITE (6,301  INO(I)  ,W(1 ,1  A) ,PHl ( I) ,W(1, 15) .COHESlI) >  W(I, 16) ,POISN(I) 
*,W(  1,17) .BULK (I) ,W( I, 13) ,S£CMOO(I> , H(1, 19) ,YOUMOO( I) ,W( 1,20)  , 

* SHE MOO (I),W(I,21),SEIVEL(I),n(I,22) , SHE VOL ( I) , 

*W(I,2A) .ATTLL(I) ,W(I,25) .ATTPI(I) ,W(I,26) , CORE ( I)  ,  W ( 1 , 27)  , 

•SMELT (I) ,W( 1,28) , VAPOR( I) , (RMKMP { I , J) , J=1 , 2) 

IF(I.EQ.NOATA)GO  TO  600 
NLINES=NLINES*2 
IF(NLINES.EQ.59)G0  TO  A00 
GO  TO  560 
600  J=0 

IF (NCARO(J*l) . EQ. 0 ) GO  TO  700 
620  WRI TE (6, AO ) 

NLINE3=2 
630  J*J*l 
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IF  (  NCARO  (  Jl  .GT.1IGO  TO  650 

WRITE  (6, 50 1  RNTE<JI  .  <  TEXT  (J,KI  ,K  =  1,19I 

NLINES=NLINES*2 

GO  TO  670 

650  WRITE  (6i60i  (TEXKJ.K)  ,K=1,1C.'> 

NLINES=NLIN£S*1 

670  lF!NCARO(JH)  .EQ.OJGO  TO  700 
IF(NLINES.GE.60)GO  TO  620 
GO  TO  630 

700  WRITE(6,70t NOATA 

10  FORMAT(»0*,2I4,6A4,2a3,  2A4 , 2X , 3 A4 , A  2 , FI  3. 2 , 3F8 . 2, F  13. 2  ,  A2 , 13,  A2  , F6 
*.2,A4) 

2  0  FORMAT (*0*.  I4,4A4,A2,F6. 0,  2A4 , 2  A3 , 1 4 ,  A2 ,  F S  . 2 ,  A2 ,  F6 . 1 ,  A  2 ,  F5 . 1 , 2  ( A2  , 
*F6.0),A2,F7.0,42,F6.4,A2,F8.0,A2,F7.0,A2,r6.4> 

30  FORHAT!*0*,I4,A2,F5.l,A2,F7.0,A2,F4.2,A2,F5.2,3(A2,F9.0) ,2(A2,F7.0 
*> ,2 (A2.F5. 1)  ,A2,F4.0,2(A2,F7.ll)  ,A4,A3) 

40  FORMAT (*1*»36X,*N  0  T  E  S*) 

50  F ORMA T(*Q*,F4.C,18A4,A3) 

60  FORMAT!*  » , 4X , 18 A4 , A3 ) 

70  FORMAT ( *1  *  «  *  THE  NUMBER  OF  OATA  ITEMS  =»,I5) 

RETURN 

ENO 


SUBROUTINE  COHEAO 

C  THIS  SUBROUTINE  PRINTS  THE  HEAOINGS  FOR  THE  CRATER  OATA  LISTING. 
WRITE!6,10! 

10  FORMAT  <*1**53X**C  RATER  OAT  A *// 1 8X. • IOENT IF ICATION* , 24X , 
••EXPLOSIVE  OATA*, 8X, ‘OEPTH* , 10X, ’APPARENT  CRATER  0 1  ME NSI ONS* , 7 X , * 

•RMK»/2X,  » . - . - . - . 

♦ - - -  0-‘  . . . . . 

*-  — -*) 

WRITE  (6,20) 

20  FORMAT!*  *,*  EVT  REF  SERIES/SHOT* , 5X , »S ITE  OATE  MEDIUM  TYPE 

*  YIELO»,5X, »EQUIV  WT  OURST  RAOIUS  DEPTH* , 5X ,* VOLUHE  SLP 

»  LIP  HI  SEEV2X  ,  *NO.  N0.*,7X,*NAME*,14X,*M0  YR*  ,  2  8  X ,  *  LBS- TNT 

*  FT  FT  FT»,7X,»CU  FT  OEG  FT  NTE»/2X,* -  -  -- 

, _ _ _ _  _ _ _ _ 

*  - - - - - - •) 

RETURN 

ENO 


SUBROUTINE  HPHEOl 

C  THIS  SUBROUTINE  PRINTS  THE  HEADINGS  FOR  THE  FIRST  HALF  OF  THE  HATERIAL 
C  PROPERTY  LISTING 
HRITEI6, 331 

30  F0RHAT(*1*,43X,*M  ATERIAL  PROPERTY  OAT  A*//,* 

♦  */10X,*SHOT  I0ENT.*,14X, 'MATERIAL*, bX,*W 

♦T-VOL  RELATIONSHIPS*, 20 X , 'STRENGTHS  ANO  STRAINS*/2X, - - 


. - . ») 

HRITEIG.'tO) 

40  FORHAT (*  *,♦  EVT  SERI ES/ SHOT *, 5X, »ELE V  HEOIUM  USDS  REF  SP  GR 

*  ORY  UWT  HOIST  TNSLE-S  TNSLE-O  UNC  OOHP  FLE  STN  DONE  COUP  CNF  PRE 
•S  FLE  STN*/2X,*N0.*,7X,*NAM£*,9X,*FT*,11X, ‘CLASS  NO.* ,8X,*LB/CUFT 

*  PR  CT  PSI  PSI*  »GX,‘PSI  IN/IN  PSI * , 7X , *3S I  IN/IN*/ 

»2X,» . - . 

*  _ _  _ _  _ _ _  _ * 

*) 

RETURN 

ENO 

SUBROUTINE  HPHE02 

C  THIS  SUBROUTINE  PRINTS  THE  HEAOINGS  FOR  THE  SECOND  HALF  OF  THE  MATERIAL 
C  PROPERTY  OATA 
HRI TE ( 6, 50 ) 

50  FORHAT  1*0*,*  */9X,‘A0(-N  STR-STN  PARAMETERS 

*  * , 1 3X  ,* MOO ULUS  VALUES', LOX, ‘HAVE  VELOC I T IES* , 2 X , *  A T I BRG  LIMITS*, 5X 

*,*  INTERNAL  ENERGIES  REMKS*/’  . 

» _ _ _ _  _ _ 

•CORE  — . . . *) 

WRITE (6, 60 1 

60  FORMAT!*  *,*  EVT  PHI  COHESION  POISN  BULK  SECANT  YOUNGS* 

* ,bX , ‘SHEAR  SEISMIC  SHEAR  LL  P I  *  ,  3X, *RECOV * , 3X , *  MELT 

*  VAPOR  SEE*/*  NO.  OEG  PSI  RAT  10  FACTOR  PSI  *,7X,* 

‘PSI*, 8X, *PSI», 7X,*FP3  FPS* ,4X, *PR  CT  PR  CT  »R  CT  MPSI/CIN  MP 

*SI/C  IN  N0TEV2X,* . - . . 

*  _ _ _ _ _ _ _ _ ......  _ _ _ _ 

. - - »1 

RETURN 

ENO 


SUBROUTINE  SURflT(NO) 

C  THIS  SUBROUTINE  DOES  A  LEAST  SQUARES  SURFACE  FIT  FOR  RAOIUS,  OEPTH,  ANO 
C  VOLUHE  OF  THE  APPARENT  CRATER 

0 IHENS IONSEC  ( 200  )  ,LX(4t),Lr(4l>,AU(200,40)  , COL  (200)  ,  XMEAN(  AO)  , 

'VECT (40) , XKOEF(40,40I  ,<<(1600), 3(40  ), 

'HOLXK(40,43) , R(40, 40), UNAT (40,43) ,OVAR(3) 

COMMON  UUM (  2  00  ,78) , INO<  2 00 ) , I  BUM ( 20 0 , 3 ) 

COMMON  /ATE/  UAN(40) ,ANS(4Q) , EXPO 
EQUIVALENCE  (XKOEF.XO 

DATA  80.CC, 00, EE , S HAl/4 HXKOF , 4HPRNT , 4 HUNI T , 4HC0RK, 1 .E-06/ 

OATA  OVAR/6HRAOIUS,bH  OEPT  H, brtVOLUHE/ 

C  M=1  FOR  RAOIUS,  2  FOR  DEPTH,  3  FOR  VOLUME 

C  L=1  FOR  EXPO=  0.250,  i  FOR  0.292,  3  FOP.  0.3125,  4  FOR  0.333,  5  FOR  VAR  EXPO. 
C  KK=1  FOR  3  PARAMETER  EQN.,  2  FOR  4  PARS.,  3  ANO  4  FOR  2  MAT.  PROPS., 

C  5  ANO  6  FOR  3  MAT.  PROPS.,  7,  8,  9,  ANO  10  FJR  4  MAT.  PROPS. 

00  900  M=l,3 
L  =  3 

00  900  KK  =  4 ,  10 

IF(KK.EQ.6.0R.«.£Q.MG0  TO  900 

CALL  EQN(COL,AU,NO.NP,L,CNORH,XMEAN,R,SEC,KK,M> 

33  HRIT£(6,1|EXP0,0VAR(H) 

1  F0RMAT(*l','£XP0NENT=',cb.4,'  FOR  ',A6,'  EVALUATION*) 

CALL  HATPRN(R,NP«£E,CC) 

CALL  COEFFICOL.AU, VECT, XKOCF.NO, HP, M) 

HRITEI6.2) 

2  FORMA  T ( ' 1 '  /  *  PRESENT  CONTENTS  OF  vECT  MATRIX*/) 

HRI TE (6,51  (COL ( 1 1 , 1  =  1 ,N°) 

5  FORMAT (*0 ' , F 12 . 4 ) 

CALL  MATPRN(XKOEF,NP,0D,CC> 

00  100  in.NP 
00  100  J=t,NP 
HOL  XK ( I , J)  =XXOEr ( I , J) 

100  CONTINUE 

CALL  MATINV(XK,LX,L1 ,NP,40,SMAL) 

CALL  MULT(U.MAT,XKOEF,HOLXK,NP) 

CALL  HATPRNIUHAT ,NP,9U, CC) 

150  CALL  SIMALT (XKOEF, VECT.ANS.NP) 

HRI TE (6, 8) 

t)  FORMA  T  ( *  1  *  /  *  NORM.  COEFS.  OF  THE  PREDICTION  EQUATION') 

CO  200  K'l.NP 
200  HRI TE (6, 13)  K.ANS(K) 

10  FORMAT ( 'O', 'N(', 12, 'I =',F19.10) 

00  90  K=1,NP 

90  UAN(K) =ANS(K)/XMEAN(<! ’CNORM 
CALL  ABPRN ( UAN , NP, KK I 

CALL  PREOIC(COL,AU,ANS, I  NO, NO , NP , M, CNORM , S EC) 

900  CONTINUE 
RETURN 
ENO 

C"' 
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SUBROUTINE  EQN1COL , AU ,NO , NP,  L.CNORH, XHtAN.R.SEC, <<,H) 

C  THIS  SUBROUTINE  OSTERNINES  THE  TYPE  REGRESSION  ANO  VARIABLES  TO  BE 
C  USED  IN  THE  LEAST  SQUARES  SURFACE  FIT  ANO  QEVELOPS  THE  AU  ARRAY 
D I  HENS  ION  COL (NO  I , AU ( NO , AO) , XME AN (A  0! , R ( 4  0 i 40) 

REAL  HEO 

DIMENSION  HEO  (2  00)  ,  T  YPE  (6)  ,  SEC  (NO) 

OATA  TYPE/iH  R, Sri  SP.JH  S,3H  M,3H  C,3H  Cri/ 

COHHON  H(200,28) .SHOT ( 200 , A) . S It E  (20  0 , 2) , DATE <20  0 , 2)  , 

*£H£D (200 . 2  I  .EXTYPE (200) , YIELO ( 200 . 3  I , TNI  HI (200 ) , 003 ( 20 U) , 
•RAOIUS(200) ,DEPTH(200), VOL (200)  ,HTL  IP (200)  , 

*RMKCO(200) . ELEVI200) , Cl  ASS (200 , 2) 

COHHON  SPGR ( 200 ) ,UWT ( 23 0) , PHOIST (20  01  , SPTEN(200) ,01  TEN (200) , 
»UCOKP(20  0) ,UCSTN(200) ,CCOHP(2QO) .CONFP (200  ) ,CCSTN(  200) ,PHl (2001  t 
»COHES(20  0) , POISN (20  0 ) ,BULK(20  0)  ,SECHOO(2UO)  ,YOUHOO(200)  , 

•SHEHOO (200 ) , SEIVEL ( 2001 ,SHEVOL(200I ,ATTLL(200), 

♦ATTPK200)  .CORE  (200)  .ENELT  (200)  .VAPOR  (200)  , RMKHP (200 , 2) 

COMMON  INO (200 ) . HREF (200) . ISL ( 200) , HREF (200) 

COHHON  /AYE/  UAN  ( <.0  )  ,  ANS  ( A 0 )  ,  EXPO 
CNORH=0.0 

GO  10(30,31,32,33.34)  ,L 

30  EXP0=l./4. 

GO  TO  40 

31  EXPO=7./?4. 

GO  TO  40 

32  EXPO=S./15. 

GO  TO  40 

33  £XPO=l./3. 

GO  TO  40 

34  EXPO=1.0 
40  CONTINUE 

00  50  J=  UNO 

43  HURT  =  (UUPMOIST(J)/tOO.)»UHT(J))/62.43 
SHUHT=HUHT*‘ (EXPO) 

SICaU0-(UWT(J)/62.43)*(  ( 1 ,0/SPGR(  J) )  »PHOI  ST  ( J)  /  10  0  • ) 

IF (SIG.LT. J . 0)  S  I G  =  0 . 0 
VHOIST=PMOIST( J)  »UWT<  J)  /6243. 

PROS=SIG»VHOIST 
VOORAT-PROS/  (l.O-PROS) 

OESAT=VHOIST/PROS 
GO  TO  (210,220,230)  ,M 
210  YaRAOIUS  (  J)  /TNTHT  (  J)  »•  (EXPO) 

GO  TO  250 

220  Y=OEPTH(J)  /TNTHT (J) •• (EXPO) 

GO  TO  250 

230  IF< VOLCJ) .LF.O.J)  V0L(J)=0,0 

YIVOL(J) »*( 1./ 3.) /TNTHT (J) **( EXPO) 

250  X-  OOB(J)  /TNTHT (J) *»(EXPO) 

IFIY.LE. 0.000) Y=  .01 
IF ( X.LE.O. JOO) X-.01 
YY= ALOG( Y) 

XX  =  ALOG(XH.O) 

COL(J)=YY 

CNORH=CNORH*-COL(J)  /NO 

GO  TO  (110  ,  120, 130,  135,  140, 140  ,  150,  150,  150 , 150) , KK 
110  AU ( J i l ) =COL ( J) 

AU  (  J,  2 )  =X 
AU(J,3) -X*X 
NP  =  3 
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CO  TO  50 

120  AU( Ji 1) -COL ! J) 

AU(J.2)=X 

AU(J,3)=X*< 

AU (J«  4) =X*  *  3 
NP  =  4 
GO  TO  50 
130  AUlJi 11=1.0 
AA=SWUHT 
B9=0ESAT 
AU( J  t  2) =44 
AU ( J i 3 ! =03 
AU(J,4)=AA*»2 
AU ( J i 5) =80  *  *  2 
AU(J,6)=A4*3B 
NMP  =6 
GO  TO  180 
135  AU<J,1>=1.3 
AU< J,2) =SHUMT 
NMP  =  2 
GO  TO  180 
140  AU (J«l)=t«J 
AA=SHUWT 
BB=0ESAT 
CC  =  SPGR<  J) 

AU ( J t  2 ) =A4 
AU ( J i 3 ) =00 
AU(J,4)=CC 
AU(J,S)=A4**2 
AU(J,6)=BB**2 
AU(J.7)=CC*»2 
AU(J,8)=A4»00 
AU( J,9)=A4*CC 
AU ( J  > 10) =8G*CC 
NHP= 1 0 
GO  TO  180 
150  AU«Jtl)=l.O 
AA=SHUHT 
OB=OESAT 
C2=SPGR<J> 

PH£=PHI( J) 7180.0*3. 141593 

IF(KX.£Q.9ICC=T4N(PHE) 

IF(K<«£Q.lO)CC=SEIVEL(J)**(l./3.) 

A  U  (  J  i  2 )  =  A  4 

AU ( J*  3) =33 

AU ( J«  4 1 =CC 

AU (J i 5) *44  *  *  2 

AU ( J i 6) =V4POR I J/ 

AU  (  J  i  7 )  =CC  *  *  2 
AU ( J i 8) =4A  *  30 
AU(J,9)=AA*CC 
AU ( J 1 1 0) =00*02 
NHP=10 
180  MS=NHP*1 
NP=NHP*3 
00  300  HH  =  MS  ?  NP 
300  AU(J.MM)  =AU<  J,MM--NMPI  *X 

IF (KK.EQ.3«0R.K<«t0.5)G0  TO  47 
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IPIKK. £0.7)00  TO  47 
HS=NP»1 
ME  =  NP*-NMP 
00  40  0MM=HS «  MF 

400  AU(  J,MM)  =  A  U  C  J,MM-NP)  »X'  *3 
NP=NPfNHP 
47  AU ( J, 1 I -COL ( J) 

50  sec ( j) =000 iji/inthtiji" texpo 

55  00  60  K=l»  NO 

COL  <K)  =COL  ( K I  /C.NORM 
60  CONTINUE 

00  70  J-l.NP 
XHEAN(J) =0.0 
DO  70  1  =  1,  NO 

XHE  AN  ( J)  =  X  NE AN  <  J I  'AU  <  I,  J)  /  NO 
70  CONTINUE 

00  80  J= 1 , N3 
00  80  1=1, NO 

AU< I , Jl =AU ( :  ,J)/ XNCAN  ( J ) 

80  CON  1  I NUE 

CALL  CORCO (AU,R,NO,NP) 

00  85  K= 1 ,  NO 
85  AU(K,l)=I.O 
XMEANI 1>=1.3 
100  RETURN 
ENO 

SUBROUTINE  C ORC 3 ( X , R , NO , NC ) 

0  HENS  ION  X(NO,40)  ,R(40,401 
00  100  1=1, NC 
00  100  J=1 , NC 
100  R ( I , Jl =0 . 0 

00  400  J= 1 , NC 
00  400  I  =  J , NC 
TOPC  =  0 • 0 
TOP  1  =  0 . 0 
I  OP J  =  0. 0 

aor i=o.a 

DOT  J  =  0 . 0 
00  800  K= l , NO 
T  0PC  =  T  OPC  K  (  K ,  I )  '  X  (  K ,  Jl 
TOPI  =  TOPIM  (K,  I) 

TOPJ=TOPJK(K,  J) 

BOTI=BOTI'X(K, I) *'2 
aor j=  bot jfx (k, ji **E 
300  CONTINUE 

R(  i  ,J)  =  «  TOPC -TOP  I/NO'  TOP  Jl  /  SQR  T  ( ( 1)  0  T  I -TOP  I/NO'  TOP  1 1  '  (  BOT  J -TOP  J/NO 
"TOPJI  I 
R 1 J , I) =R< I , Jl 
400  CONTINUE 
RETURN 
ENO 
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S UO ROU TI N£  ft B P R H ( J A N , NP i  KK ) 

DIMENSION  UAN(40) ,8(401  ,  04(2) 

WRITE  (6, 8) 

8  FORHA T (* 1*/*  COEFFICIENTS  OF  THE  PREDICTION  EQUATION*) 
00  200  K=1,NP 
200  HRITE(6.10)  K.UAN(K) 

10  FORHAT (*  *,*C<*, 12,*) =*,F19.1Q) 

GO  TO  (310, 320,800,800, 800, SuO, 800, 600,800,800) ,KK 
310  0(1)=  EXP (U4N( 1) -  UAN(>) ♦ ' 2/ ( A . 0* Oft N( 3 ) )) 

0(2)  -UAN(J) 

0(3) =UAN(2> / ( 2*OAN( 3 ) ) 

HRI TE (6.8S) 

88  FORMAT ( *0*^ »  COEFFICIENTS  OF  THE  PREDICTION  EQUATION*) 
DO  101  MM=1,NP 
101  HRI TE (6, 1 1 >  HH.OIHH) 

11  F  ORM  A  T  ( *  0*  ,  *  0  (  *  ,  12,  *)  =•  ,Fl9.  10) 

GO  TO  800 

320  0(2)  =UAN<4) 

0T=UAN(3)/UAN(4) 

0P=  AOS(DT) 

UNR AO=BP” 2  ♦  9*  ( 3  ♦  O * UA N(  2 )  7UAN <  4 ) ) 

IF (UNRAO.GE. 0.0 1  GO  TO  755 
HRI TE (6, 709)  UNRAO 
200  FORHAT(»0*,»UNRAO=»,F2G.IO) 

UNR A0=0 . 0 

755  3SQ  T  =  SORT (UNRAO) 

04(i)=(or*asar)/3.o 
04(2) =(BT-DSQT)/3.0 
00  900  J J=l , 2 
813  0<4>=04<JJ1 
812  0(3) =OT-(2. 0*0(41 ) 

0NEG  =  AOS ( 3 (  4) ) 

0(1)=  EXPlUftNI 1) -  8(21*0(3) *UNEG* * i ) 

HRI TE (6, 55 ) 

55  F  ORMA  I  ( *  0  * / *  COEFFICIENTS  OF  TmE  PREDICTION  EQUATION’ ) 
00  191  MMM  =  1 , NP 
191  HRI TE (6,2?)  MMH, O(MHM) 

22  FORMAT  (*0*.  ’()(’,  12,*)  =*  .F19.10) 

900  CONTINUE 
800  RETURN 
END 


SUJROUT INE  COEFF (COL. ftU, VECT , XKOEF, NO, NP.M) 

c  this  suoroutine  generates  ano  assembles  the  vect  and  xkdsf  matrices 

DIMENSION  COL (NO) . AU( NO ,40) , VECT (40) .XKOEF (40,49) 

ADA  T  A  =  NO 
00  10  1  =  1, NP 
VECT ( I)  =0. 

10  CONTINUE 

00  30  1=1, NP 
00  25  J= l , NO 

VEC  T  ( I)  =VEC  T  ( I )  *COL  (  J  )  *  A U(  J,  I  > 

25  CONTINUE 

VEC  T(I) =VEC T ( I ) / AOATA 
30  CONTINUE 

00  20  J  =  1 «  NP 
00  20  K=1,H? 

XKOEF  U.K) =0  .0 
20  CONTINUE 
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00  40  1  =  1, NP 
00  40  J=1,NP 
00  30  K= 1 , NO 

XKOEF  (I, J)  =AU(<,  I)  »AJ(K,o)  ’XKOEF (J,J> 
35  CONTINUE 

XKOEF ( I, J) =XKOEF(I,J)/AOATA 
40  CONTINUE 
50  RETURN 
ENO 


SUBROUTINE  PREOICICOL ,A U , ANS , I E V , NO , NP , H, CNORH, SEC l 
C  THIS  SUBROUTINE  CALCULATES  THE  RAOIUS,  0E°  T  H ,  OR  VOLUME  USING  THE 
C  EHPIRICAL  EQUATION  GENERATED  ANO  COMPARES  THESE  VALUES  WITH  THE  ACTUAL. 
0 1  HE NS I ON  lEV(NO) ,COl (HO) ,AUtNO,40)  ,ANS tNP) ,0VAR(3)  , SEC (NO) , PARTE) 
OATA  0VAR76HRA0IUS, bH  OEPTH, bHVOLUME/ 

NRlTEIb,  10)  OVAR(M) 

lti  FORMAT ( *  1  *  ,  4  b ,  *  PREOICTIOH  ANO  EVALUATION*) 

HRI TET6.20) 

20  FORHAT(»-*,10X,*PREOICTEO  VALUE *, 12X ,» AC TUAL  VALUE* , 12X, 'PERCENT  c 
»RROR*,12X,*RESIOUAL*.  9X, ’EVENT  NO.*,bX,*SC  000*) 

Ch£  4N=0 . 0 
00  90  ' 1 , N 0 

90  CMt  AN=CH£  AN  f  EXP  (COH  I)  *CNORH) 

RTOP=0. 0 
R30T  =  0.0 
NHNES=0 
00  200  1=1, NO 
PREVAL  =  3.G 
00  100  J= 1 , NP 

PR£VAL=PREVAl*AU(I . J) »ANS( J) 

100  CONTINUE 

PREVAL=PREVAL’C.NORM 
ACTVAL=COL (I) *CNORM 
PREVAL=  EXP(PP.EVAL) 

ACTVAL=  EXP(ACTVAL) 

IF ( ACTVAL.EQ.O.O) AC T V AL  =  1 . OE -  3 
RESIO=PREVAl-ACTVAL 
error=resio/actval*  ioo. g 

HRIT£(6,33IPREVALiACTVAL,ERROR,RESIO,IEV(I)  ,  SECII) 

30  FORMAT (»  *,6X,FlV.6,3X,Fl7.6,15X,F8.2,6X,F17.b,9X,I5,2x,Fl5.4) 

RT 0P=RT0-> )  (PREVAL-CNORN)  **2 
R80T  =  R30TKACTVAL-CN0RA)  ”2 
NL I NES=NL INES  *  1 
IF(NLIHES.LT.50)GU  TO  200 
NLIHES=0 
HRI TE (b, 50 ) 

50  FORMAT ( “ t* ) 

200  CONTINUE 

RSQ=RTOP'RBOT*100. 

S=SQRT (ROOT '  (NO-1) ) 

HRI IE (b, 40  I  RSQ.OVAR (ri) • 3 

40  FORMATv»0*,*MULT  CORR  COEF  = • , F9 • 2/ *  STAND.  DEVIATION  FOR  *,Ab,*=* 

’ ,  F 1 2 . 4  < 

RETURN 

ENO 
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SUBROUTINE  MAUNV  ( A  ,  IRON,  TOOL  ,  NP,  NOI M ,  SML  ST) 
DIMENSION  A  (  1600)  ,  IROHI  4  1)  ,  IC0L141) 

N=NP 


NPl=N'l 

GUNA 

136 

00  5  1=1, N 

GUNA 

137 

ICOLI  I)  =1 

GUNA 

138 

5 

IROH ( I)  =  I 

GUNA 

139 

00  75  ITER=1,N 

GUNA 

140 

MAXR= ITER 

GUNA 

141 

HAXC=1 

TE.1P=  AOS  (A  l  HAXR)  ) 

GUNA 

142 

L IHI TC=NP1 - 1  TER 

GUNA 

144 

00  15  I=IT£R,N 

GUNA 

145 

00  15  J=l, LIMITS 

GUNA 

146 

IJ=(J-1)  'NOIMH 

GUNA 

14  7 

IF  I  TEMP- (  AQSCA (I J) ) ) i  10,10,15 

GUNA 

148 

10 

HA  X  R= I 

GUNA 

149 

MAXC=J 

GUNA 

150 

TEMP=  AQS ( A ( IJ) I 

LUNA 

151 

10 

CONTINUE 

GUNA 

152 

SML  ST =-0.0 

GUNA 

153 

IF  (TE.'P-SNLST)  20,23,25 

GUNA 

154 

20 

IRON (NP1 ) =  I  T  ER 

GUNA 

155 

200 

KRI TE (6, 200) 

FORMAT ( * 0 * , ’TrilS  IS  A  SINGULAR  MATRIX  AND  IT  RILL  NOT  INVCRT*) 

oUNA 

156 

25 

IF  (MAXR-ITER)  30,40,33 

GUNA 

159 

30 

00  35  J=1,N 

GUNA 

160 

HA  XRJ=  < J- 1 ) 'NOIM'MAXR 

GUNA 

161 

ITJ=(J-1)  'NOIMHTER 

GUNA 

162 

IEMP=A(MA<RJ) 

GUNA 

163 

A(MAXRJ)=A(ITJ) 

GUNA 

164 

35 

A(ITJ)=TEMP 

GUNA 

165 

ITcMP=IROR (MAXR) 

GUNA 

166 

IROH(MAXR)  =  I  ROR  (  HER) 

GUNA 

167 

IROH<  I  TE.RI  =  I  TEMP 

GUNA 

168 

40 

IF  (MAXC-t)  45,55,45 

GUNA 

169 

45 

00  50  1=1, N 

GUNA 

170 

IMAXC  =  (MAXC- 11  'NOI M* I 

GUNA 

171 

TEMP  =  A(  I) 

GUNA 

172 

A  (  I )  =  A  ( I MA 

GUNA 

173 

50 

A( IMAXC) =TE5P 

GUNA 

174 

ITtMP=lCOL (MAXCI 

GUNA 

175 

ICOL(MAXC) s I  COL ( 1 1 

GUNA 

176 

ICOL  ( 1)  =  I  TEMP 

GUNA 

177 

55 

T£MP=A(ITER) 

GUNA 

178 

IT6MP=IC0L (1) 

GUNA 

1  79 

00  60  J=?,N 

GUNA 

180 

IT JMlr(J-r  I  'NOIM'ITtR 

GUNA 

181 

IT  J=  (J-l)  *NOIMU  TER 

GUNA 

182 

A  (  I  TJM1)  =A  » I  TJI  /TEMP 

GUNA 

183 

60 

ICOUJ-l)  =  ICOL(J) 

GUNA 

184 

I TN  =  <  M-l » ’NOIM' I  TER 

GUNA 

185 

A(  ITN)=1. OFTEN® 

GUNA 

186 

ICOL(N)  =IIEHP 

GUNA 

187 

00  75  1=1, H 

GUNA 

188 

IF  <  I  - 1 T !.  R 1  65,75,65 

GUNA 

189 

65 

TEMP-A ! I 1 

GUNA 

190 

00  70  J=  2 , N 

GUNA 

191 
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IJM1*(J-2)»N0IM*I 

GUNA  192 

IJ*IJ-1J*N0IM*I 

GUNA  193 

irjHi»u-2i  ’noikmter 

GUNA  194 

A(lJMi)«A(IJ)-A(tTJHll*TEHP 

GUNA  195 

70 

CONTINUE 

GUNA  196 

INMN-1>*M«IM*I 

GUNA  197 

ITN*(H~1)*N0IN*ITER 

GUNA  198 

A(IN)»-(TEMP*A(ITN)  1 

GUNA  199 

75 

CONTINUE 

GUNA  200 

00  100  I *  1 1 N 

GUNA  201 

DO  50  J-’rn 

GUNA  202 

IF  <1 ROMIJI-Il  80,85,80 

GUNA  203 

80 

CONTINUE 

GUNA  204 

85 

IF  (I-J)  90,100,90 

GUNA  205 

90 

OO  95  L*1,N 

GUNA  206 

LI*  (I-1)*N0IM*L 

GUNA  207 

LJ*<J-1)  »N0IM*L 

GUNA  208 

TEMP*A(LI) 

GUNA  209 

AIL  I) *A(LJ> 

GUNA  210 

95 

A(LJ)«TEM» 

GUNA  211 

IRON ( J) *  IRON ( I ) 

GUNA  212 

100 

CONTINUE 

GUNA  213 

00  125  1*1, H 

GUNA  214 

00  105  J*I,N 

GUNA  215 

IF  <IC0UJ>-I>  100,110,105 

GUNA  216 

105 

CONTINUE 

GUNA  21? 

110 

IF  II-J)  115,125,115 

GUNA  218 

115 

00  120  L*l,N 

GUNA  219 

lL*IL-i)*N0IMM 

GUNA  220 

JL*  (L-l) *N0IM*J 

GUNA  221 

T£NP*A  III) 

GUNA  222 

A 1 1 L ) *A ( JL 1 

GUNA  223 

120 

AIJD-TEMP 

GUNA  224 

ICOL  ( J)  *IC0L  ( I ) 

GUNA  225 

125 

CONTINUE 

GUNA  226 

IROHINPl)  *0 

GUNA  227 

RETURN 

GUNA  228 

ENO 

GUNA  225 

SUOROUTINC  MATPRN  I  A  ,  NP ,  8  ,  RITE  ) 

OIHENSION  A  l*»0,**0> 

N«NP 

5 

FORMAT  (»0*//) 

10 

FORMAT<*0*,11IE12.3)) 

22 

FORMAT  |*l*  ,  2 4H  PRE5EHT  CONTENTS  OF 

»  A  * 

SUP  1277 

« 

»  12H  MATRIX  ,  t  ) 

OATA  YA2/AHPRNT/ 

IF  1  RITE  .NE.  TA Z  l  GO  TO  £00 

SUP  1280 

MRITE  (6,22)  8 

SUP  1285 

IF  I  N.CT.ll  )  GO  TO  6) 

SUP  1286 

00  50  I«1,NP 

MRITE  (6,10)  I A I  I ,  J)  ,  J*  1 ,  NP) 

50 

CONTINUE 

SUP  1289 

GO  TO  500 

SUP  1290 

60 

Ml  «  N/11M1-10 
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00  200  K  =  1  ,  N1  ,  11 
«  :  <  I  II 
00  100  I  =  1  ,  N 

WRITE  (6 i 1 0 )  (  4(1, J)  ,  J=<  ,  «) 
100  CONTINUE 
WRITE  (6,5) 

200  CONTINUE 
<=KU1 

00  000  I  =  1  ,  N 

WRITE  (6,13)  (  A(I,JI  ,  J=<  ,  N  i 
300  CONTINUE 
600  CONTINUE 
50  0  RETURN 
ENO 


SUOROUTINE  SIMALI  (A.O.C.NP) 
OIHENSIOM  4(40,40),  0 (♦ 0 ) , C ( 4 0 ) 
M  =  MP 
N-NP 

00  200  I  s  1  ,  H 
C(I)  =  0.0 

200  CONTINUE 

00  500  I  *  1  ,  M 
00  400  J  *  1  ,  N 

C(  1)  =  C  ( I)  *■  A  (  I  ,  J)  »  0(  If 

40 u  CONTINUE 
500  CONTINUE 
RE  TURN 
ENO 


SUP  1292 
SUP  1293 
SUP  1294 
SUP  1295 
SUP  1296 
SUP  1297 
SUP  1299 

SUP  1300 
SUP  1301 
SUP  1302 
SUP  1304 

SUP  1306 


E.UNA  124 
uUNA  125 
GUNA  126 
GUNA  127 
OUNA  128 
uUNA  129 
GUNA  130 
GUNA  131 
GUNA  132 
GUNA  133 


SUJROUTINE  HULT(A,3,C,NP) 

OIHENSIOM  A(4Q, 401, 0(40, 401, 0(40,431 


K*NP 
M  =  NP 
N  =  NP 

00  100  1=1, K  OUAR  723 

00  200  J  =  1 , N  OUAR  724 

>3  =  9.0  OUAR  725 

00  300  1*1. H  QUAR  726 

Q  =  0  *■  9(1, U  *C (L, J)  QUAR  727 

300  CONTINUE  OUAR  728 

A ( I , J)  =  Q  OUAR  729 

200  CONTINUE  QUAR  730 

100  CONTINUE  QUAR  731 

RETURN 
ENO 
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FIG.  10.  CRATER  VOLUME  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  DESERT  ALLUVIUM 


120 


□EPTH/W**5/16  RRDIL6/N*»5/16 


FIG.  n.  CRATER  RADIUS  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  SAND 


FIG.  12.  CRATER  DEPTH  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  SAND 
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(VQL**i/3J/W**5/l6 
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FIG.  23.  CRATER  RADIUS  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  PLAYA  (AIR  VENT  SERIES) 
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FIG.  24.  CRATER  DEPTH  AS  A  FUNCTION  OF  CHARGE, 
DEPTH  FOR  PLAYA  (AIR  VENT  SERIES) 
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RRDIUS/N*»5/16  IV0L*«1/3)/Wk*5/16 


FIG.  25.  CRATER  VOLUME  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  PLAYA  (AIR  VENT  SERIES) 


FIG.  26.  CRATER  RADIUS  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  PLAYA  (TOBOGGAN  SERIES) 
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IV0L**1/3J/W**5/I6  0EPTH/W**5/ 16 


FIG.  27.  CRATER  DEPTH  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  PLAYA  (TOBOGGAN  SERIES) 


FIG.  28.  CRATER  VOLUME  AS  A  FUNCTION  OF  CHARGE 
DEPTH  FOR  PLAYA  (TOBOGGAN  SERIES) 
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APPENDIX  IV 

GENERAL  EQUATION  COEFFICIENTS 


TABLE  4.  GENERAL  EQUATION  COEFFICIENTS  FOR  RADIUS 


Using  S  and  Gc 


Using  y^6,  S,  G$  and  Ey 


C(  6)  = 

C(  7)  = 

Cl  8)  = 

C(  9)  = 

C(10)  = 

Clil)  = 

C  ( 12)  = 

C  c 13)  = 

C  ( 14)  = 

C ( 15) = 

Cl  16)  = 
C<17)  = 
Cll8)  = 

C(  19)  = 
C(20)= 
C(21)  = 

CT 22) = 

C  ( 23)  = 

C  l  24)  - 
C(25)= 

C (26) = 

C ( 27 ) = 
C128)= 

C  (29)  = 
C(30)  = 

4  ■  95-9 

63%  within  ± 
88%  within  ± 


-31.8910003853 
53.7596378609 
-8.9711680343 
.5549490163 
-30.5239196038 
-3.373 3199735 
-.9415056763 
14.0810150793 
5.0934471223 
-1.9052946210 
-61.0822676775 
10.5127003945 
-2.7082997779 
41.4257955423  i 
44.6352964490 
4.3584334973 
1.7054916517 
-9.6728410805 
-42.7140631702 
3.8425311344 
71.0614236931 
-.8820397422 
1.7405271634 
-53.1325358655 
-30.3299093539 
-1.4850841324 
4.1108960133 
5.4834047720 
26.7268983051 
-2.6092817849 


Cl  1)  ? 

Cl  2)  = 

C(  3)- 
C(  4)  = 

Cl  5)  = 

C  (  6)  - 
C  (  ?>  = 

C(  8)  = 

Cl  9)  = 

Cl  10)  = 

c  c  ID  = 

C ( 12) - 
C  ( 13)  = 
C(14>= 

C  ?  15)  - 
C  ( 16)  = 

C ( 17) = 

C  ( 1 8)  = 

C  ( 19)  = 

C  (  20)  = 

C  (  2 1)  = 

C  ( 22)  = 
C(23)  = 

C  C  24)  = 

C  ( 25)  = 

C ( 26) = 

C  ( 27)  = 

C  { 28)  = 

C (29)= 

C  (  3  0 )  = 

C  (31)  = 

C  (  32 )  = 

C ( 33) = 
C(34)= 

C  ( 3  5 )  = 

C  ( 36)  = 

C  l  37)  = 

C  (  38)  = 

C  ( 39)  = 

C  ( 4  0)  = 

R?  =  97  • 7 

72%  within 
91%  within 


4.7431667917 
-312.2659209015 
85 . 1020403634 
130.7525251740 
2.2953850824 
-.0000081432 
-52.3758442797 
-12.4911945922 
122.363631315? 
-27.8147254654 
-141.9455105232 
323.0631524575 
-196.6141937964 
-251.5265120343 
-153.2196787513 
.0000127300 
96.2318031404 
126.6088317324 
-199.3081650433 
18.4750255695 
112.7398253179 
-551.6962384152 
124.3090113793 
156.6435857970 
184.2976000120 
-.0000075941 
-48.5622042860 
-135.4375636625 
65.7358828227 
13.9867726703 
-13.4793572557 
112.2478478970 
-24.0214629682 
-38.9836998015 
-59.5632559447 
.0000004885 
7.3256515643 
39.1585126325 
5.2268297684 
-8.6972701637 


±  10% 
=  20% 


J 


TABLE  4.  GENERAL  EQUATION  COEFFICIENTS  FOR  RADIUS  (CONTINUED) 


Using 

,  S,  tan  4>  and  Ev 

C  <  i)  = 

18.8992452280 

C(  2)  = 

-31.8938433251 

C(  3)  = 

20.3800253449 

C(  4)  = 

-8.398280971b 

C(  5)  = 

13,2362688240 

C(  6)  = 

-.0000027090 

C  (  7)  = 

.4473329001 

C(  8)  = 

-17.5521678930 

C(  9)  = 

6.8566842352 

C(10)= 

1.5946340224 

C(ll)  = 

-199.0532997153 

C ( 12) = 

323.2985371981 

C  ( 13)  = 

-116.5369735224 

C  ( 14)  = 

64.5695519242 

C(15)  = 

-128.5006629243 

C  ( 16)  = 

-.0000105683 

C  ( 17)  = 

-4.069150625b 

C  (18)  = 

94.5027765299 

C  ( 19)  = 

-51.7359125275 

C  ( 20)  = 

-1.0595754099 

C  ( 21)  = 

225.6633484300 

C  ( 22)  = 

-361.8623322793 

C  ( 23)  = 

125.7080514573 

C  (  24)  = 

-76.6928891973 

C  ( 25)  = 

140.4383152831 

C  (26)  — 

.0000176320 

C  (27)  = 

2.5802223123 

C  ( 28)  = 

-96.9396721804 

C(29)  = 

66.0886604129 

C  (  30 )  = 

-6.9956496680 

C  ( 31)  = 

-68.9475372933 

C  (32)  = 

109.5418270782 

C  ( 33)  = 

-37.8928184431 

C  ( 34)  = 

24.9124062737 

C(35)  = 

-41.8140082327 

C  ( 36)  = 

-.0000078407 

C ( 37) = 

-.3086871023 

C  (38)  = 

28.281851447b 

C ( 39) = 

-22.3678557147 

C  (  40 )  = 

3.8323586654 

=  98.1 
m 

76%  within 

±  10% 

94%  within 

±  20% 

Using  y5/16>  S,  c1/3  and  Ey 


C(  1)  = 

29.3980672435 

C(  2>  = 

17.2560192731 

C  (  3)  = 

10.020396317b 

C(  4)  = 

-4.4695237570 

C(  5)  = 

-33.2041384178 

C(  6}  = 

-.0000024075 

C(  7)  = 

-.0027363973 

C  (  8)  = 

-6.6328662697 

C(  9)  = 

3.6283822224 

C(10)  = 

-.1239487349 

C  ( 11)  = 

-122.3373482461 

C ( 12) = 

-7.4095270897 

C  ( 13)  = 

-46.8808333927 

C  1 14)  = 

14.6150371805 

C ( 15) = 

89.2867812633 

C ( 16) = 

-.0000056701 

C ( 17) = 

.0232794163 

C  ( 18)  = 

38.5678057694 

C  (  19)  = 

-12.2559501813 

C  ( 20)  = 

.0819726642 

C  ( 21)  = 

95.1305302502 

C  ( 22)  = 

38.7459131353 

C ( 23) = 

22.0664331673 

C  ( 24)  = 

-13.128434647b 

C  ( 25)  = 

-96.2455129446 

C  ( 26)  = 

.0000127685 

C  ( 2  7)  = 

-.0114900349 

C  ( 28)  = 

-13.6140033792 

C  ( 29)  = 

10.8696441627 

C  ( 30)  = 

-.3975240262 

C  (  3 1)  = 

-24.6729413515 

C  (  32)  = 

-13.6655896113 

C  (33)  = 

-.8950301232 

C ( 34) = 

3.5235047360 

C  (  3  5)  = 

27.6217414053 

C  ( 36)  = 

- .0000066210 

C  (37)  = 

-.0011715210 

C  (  38)  = 

-2.1689813563 

C(39)  = 

-2.8343303863 

C  (  40)  = 

.2215007433 

RnZ,  =  97.7 

73%  within 

±  10% 

95%  within 

±  20% 

'3* 


TABLE  5.  GENERAL  EQUATION  COEFFICIENTS  FOR  DEPTH 


Using  y 

5/16,  S  and  G$ 

Using  y5//16,  S,  Gg  and  Ey 

C<  1)  = 

53.3958631782 

C  (  1)  = 

-43.6955704922 

C(  2)  = 

-127.9226404562 

C  (  2)  = 

-263.0260424584 

C<  3)  = 

20.8301124395 

C(  3)  = 

76.0756460918 

C(  4)  = 

17.9316706900 

C(  4)  = 

143.4022465293 

C<  5)  = 

-48.4457154475 

C(  5)  = 

14.6854054001 

C  <  6)  = 

-.5180822422 

C  (  6)  = 

-.0000018621 

C  (  7)  = 

-24.1731506505 

C  (  7)  = 

-47.1730917157 

C<  8)  = 

15.8709481683 

C  (  8)  = 

"15.6964256283 

Ci  9)  = 

91.6107413501 

C(  9)  = 

90.9027550405 

C<10)= 

-15.2484978767 

C  ( 1 0)  = 

-22.2985600530 

C(il)  = 

-449.1867152981 

C(lt>  = 

87.0048968184 

C(12)= 

368.5656181405 

C ( 12) = 

660.3121352892 

C  ( 13)  = 

-44.6720935833 

C  ( 13)  ~ 

-126.3905248215 

C  ( 14)  = 

166.2582151665 

C ( 14) = 

-354.1566836962 

C  ( 15)  = 

66.2414675220 

C  ( 15)  = 

-82.2950451513 

C  ( 16)  = 

-3.7716040549 

C ( 16) = 

.0000200377 

C  ( 17)  = 

17.9221947624 

C ( 17) = 

112.9332700877 

C  (18)  = 

11.5629558875 

r  u  «i  = 

90 -8534757507 

C  ( 19)  = 

-205.3794622239 

C  ( 19)  = 

-193.3862745513 

C  C  20)  = 

13.0631609205 

C  (  2 0)  = 

8.2417379757 

C  ( 21)  = 

254.9573454245 

C (21) = 

-598.9320341567 

C  ( 22)  = 

-120.7791721681 

C  ( 22)  = 

-523.9584928841 

C ( 23) = 

14.2414607555 

C ( 23) = 

41.3421692583 

C  (24)  = 

-134.8771751525 

C  ( 24)  = 

689.3555775099 

C ( 25) = 

-51.6017664322 

C (25)= 

107.4665743070 

C ( 26) = 

1.9569970889 

C ( 26) = 

-.0000353683 

C ( 27) - 

2.9310692111, 

C  ( 27)  = 

-160.9998103999 

C  ( 28)  = 

-.1726515931 

C  ( 28)  = 

-93.9073485683 

C ( 29) = 

94.4689354767 

C ( 29) = 

118.7384287103 

C ( 30) = 

-6.0438470627 

C  (  3  0 )  = 

25.6274231773 

R2  =  93  6 

C  ( 31)  = 

342.2023179307 

m 

C ( 32) - 

175.6525289593 

39%  within 

±  10% 

C (33)= 

.3384907975 

65%  within 

±  20% 

C  (  34)  = 

-340.1575068641 

C  ( 35)  = 

-53.7051347245 

C  ( 36)  = 

.0000144483 

C ( 37) = 

71.4063952381 

C (38)= 

33.5127395672 

C  ( 39)  = 

-23.9273648252 

C ( 40) = 

-15.1138557381 

R2  =  96.0 
m 

47%  within 

±  10% 

72%  within 

±  20% 
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TABLE  5.  GENERAL  EQUATION  COEFFICIENTS  FOR  DEPTH  (CONTINUED) 


Using  y5/ 16 

,  S,  tan  $  and  Ey 

Using  S,  and  Ey 

C  (  i)  = 

17.1253607120 • 

C(  1)  = 

21.7086303701 

C<  2)  = 

-37.3542155845 

C  (  2)  = 

-42.3323530877 

C(  3)  = 

12.5105364463 

C  (  3)  = 

14.6236850663 

C(  4) - 

8.5029333445 

C  (  4)  = 

.1510212049 

C(  5)  = 

14.7403694963 

C  (  5)  = 

22.8611299509 

C(  6)  = 

.0000029697 

C(  6)  = 

.0000028219 

C(  7)  a 

-7.2058081245 

C  (  7)  = 

.0096628854 

C(  8)  = 

-8.4253426900 

C(  8)  - 

-14.1319628653 

C(  9)  = 

3  .3999072147 

C  (  9)  = 

-.4996236569 

C (10) = 

-1.8629595715 

C (10)= 

.2047476496 

CM1)  = 

-193.3918309651 

C(ll>= 

-85.5758975799 

C ( 12)  = 

330.9456858185 

C ( 12) = 

112.7124014239 

C  (13)  - 

-97.2329269833 

C  ( 13 )  = 

-25.8772673551 

C  ( 14)  = 

33.0550538845 

C ( 14) = 

2.6748938293 

C ( 15) = 

-134.3558827033 

C ( 15) = 

-45.0963840333 

C ( 16) = 

-.0000253473 

C  ( 16)  = 

-.0000128230 

C  ( 17)  = 

4.4905080385 

C ( 17) = 

-.0277178934 

C (18) = 

72.4656077635 

C  ( 18)  = 

32.4942151045 

C ( 19) “ 

-37.2785637467 

C  (19)  = 

-.9376620835 

C (20) = 

7.6485959924 

C  ( 20)  = 

-.7761191183 

C  ( 21)  = 

219.4717738692 

C  ( 21)  = 

90.0293948233 

C ( 22) = 

-361.0416787948 

C ( 22) = 

-86.5807061313 

C ( 23) = 

113.5789573079 

C ( 23) = 

-16.6599345535 

C ( 24) - 

-63.5462262740 

C (24)= 

-4.1166825850 

C  ( 25)  = 

144.1389355682 

C ( 25) = 

27.0640857132 

C(26)= 

.0000298469 

C (26) = 

.0000219664 

C(2  7)  = 

2.6975157247 

G  ( 27)  = 

.0528352719 

C ( 28) = 

-83.6172087730 

C  ( 28)  = 

6.2005227223 

C  (  29)  a 

54.6675138674 

C ( 29) = 

1.3423129962 

C  ( 3  0 )  = 

-12.3071777703 

C  ( 30)  = 

.4649869039 

C ( 31 ) = 

-72.7217500139 

C  (  3 1)  = 

-31.6097537754 

C ( 32) = 

116.2322451501 

C  ( 32)  = 

32.5439752740 

C ( 33 ) = 

-37.4935916521 

C  (  33)  = 

11.4732120724 

C ( 34) = 

25.5528645197 

C  (34)  = 

1.1107085605 

C ( 35) = 

-45.3802628657 

C  (35)  = 

-12.4902317353 

C  ( 36)  = 

-.0000108194 

C  (  3  6)  = 

-.0000104873 

C  (  37)  - 

-1.3321194171 

C  (  37)  = 

-.0264307464 

C(38)= 

27.3519570403 

C  (  3  8  >  = 

-8.2155628902 

C ( 39) = 

-21.3464789174 

[  C  (  39)  - 

-.0050793512 

C  ( 40)  = 

4.9772599019 

!  C(40)  = 

-.0721852584 

9.1  =  94.5 
m 

43%  within 

±  10% 

Rf  =  97.8 

45%  within  ± 

10% 

66%  within 

±  20% 

72%  within  ± 

20% 
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TABLE  6.  GENERAL  EQUATION  COEFFICIENTS  FOR  VOLUME 


Using 

y5/16,  S  and  Gs 

Using  y5/16,  S,  G$  and  Ey 

C  {  1)  = 

-2,1752357203 

C(  1)  = 

-11.8207236195 

C<  2)  = 

13.1141215646 

C  (  2)  = 

-347.0492941075 

C  (  3)  = 

.5997600122 

C(  3)  = 

97.9246661795 

C<  4)  = 

-3.6647026294 

C  (  4)  = 

156.8518890100 

Ct  5)  = 

-30.4405669790 

C(  5 )  = 

1U. 5575210571 

C  (  6 )  = 

-3.4757992071 

C(  6)  = 

-.0000070363 

C  (  7)  = 

-3 .8080373705 

C(  7 )  = 

-58.4668518002 

c<  8>  = 

14.9439527393 

C(  8)  = 

-16.1837586229 

C  (  9)  = 

20.8239905853 

C<  9)  = 

128.5893232243 

C ( 10) = 

-5.9957649451 

C(10)= 

-30.9888871631 

C(ll)= 

-215.7667550905 

C(ii)= 

-59.3569012349 

C ( 12) = 

93.1623369373 

C ( 12) = 

875.9625480863 

C  < 13)  = 

-16.5  0  U  8692970 

0(13)= 

-209.5974746947 

C  ( 14)  = 

118.9135549093 

C  ( 14)  = 

-334,6553293740 

C ( 15) = 

33.8867065685 

0(15)= 

-147.6046367155 

C ( 16) = 

2.7663678702 

C  ( 16)  = 

.0000132385 

C  <  1 7 )  — 

-6.3022129785 

C  ( 17)  = 

117.2272444493 

C ( 18) = 

2.7930438894 

C  ( 18)  = 

131.0815866705 

C<19)  = 

-67.2103550710  j 

C ( 19) = 

-224.9307823733 

C ( 20) = 

4.0009580199 

C ( 20) = 

21.5972214015 

C(21)  = 

146.3688958239 

C ( 21 ) = 

-134.8691756255 

C ( 22) = 

-25.5335254205 

C ( 22) = 

-617.5477326839 

C (23) = 

6.3873929627 

C ( 23) = 

122.2498533469 

C  ( 24)  = 

-98.0675342203 

C  ( 24)  = 

372.8148507687 

C ( 25) = 

-30.2435268423 

C  ( 25)  = 

163.4194029422 

C ( 26) = 

-.6234156121 

C ( 26) = 

-.0000111033 

C ( 27) = 

9.8235554535 

C ( 27) = 

-99.1900726417 

C ( 28) = 

1.2453679293 

C ( 28) = 

-133.4532348233 

C ( 29) = 

37.11d6765414 

C  (  2  9 )  = 

109 .7545480032 

C  (  3  0)  = 

-2 . 7689859994 

C (30)= 

13.5108852373 

? 

C  ( 31)  = 

109.1642310421 

R~  »  94.7 
m 

C (32)= 

152.1133991992 

59%  within 

±  10% 

C (33)= 

-21.4879680004 

84%  within 

±  20% 

C ( 34) = 

C  (  3  5)  = 

C ( 36) = 

C (37)= 

C  ( 38)  = 

C  ( 39)  = 

C  ( 40)  = 

=  97-2 

65%  within 
89%  within 

-150.2371189017 

-55.9322230211 

.0000026885 

32.7649358627 

39.6790168769 

-13.5594018533 

-9.8179269609 

±  10% 

±  20% 
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TABLE  6.  GENERAL  EQUATION  COEFFICIENTS  FOR  VOLUME  (CONTINUED) 


Using  y5/16 

,  S,  tan  and  Ev 

Using  y3^® 

,  S,  c1'3,  and  Ey 

C(  1>  = 

22.3706247625 

C<  i)  = 

20.4178474204 

C(  2)  = 

-38.9723779360 

C<  2)  = 

6.7210259119 

C  (  3)  = 

18.9610934544 

C(  3 )  = 

11 997  6882110 

C<  4)  = 

-5.7817278439 

C<  4)  = 

-2.8361992670 

C(  5)  = 

15.6413196900 

C(  5)  = 

-17.5564771871 

C(  6/  = 

-.0000019613 

C(  6)  = 

-.0000012370 

C(  7)- 

-1.4903246793 

C(  7)  = 

.0043463631 

C(  8)  = 

-15.6027125593 

Ci  8>  = 

-9.5472334444 

C(  9)  = 

7.4617827703 

C(  9)  = 

2.1150419147 

C(iO)= 

.2976374547 

C(10)= 

-.0237115989 

C(ll)= 

-205. 6266 0425Q0 

C(li)  = 

-83.7714888573 

C  ( 12)  = 

332.7916277065 

C  C 12)  = 

-1.0129467374 

C  ( 13)  = 

-107.5617213000 

C  (13)  = 

-40 . 4930003277 

C  ( 1 4 )  = 

65.3451232219 

C ( 14) = 

10.4469006002 

C ( 15) = 

-131.0986024125 

C ( 15) = 

56.2397519723 

C  ( 16)  = 

-.0000133719 

C( 16)= 

-.0000128353 

C  ( 17)  = 

-2.7122525642 

C  ( 17)  = 

-.0069710347 

C ( 18) = 

85.1580280390 

C ( 18)= 

36.2108467755 

C  ( 19)  = 

-54.5527834882 

C  ( 19)  = 

-8.0426022633 

C  ( 20)  = 

2.9801404833 

C  (  2  0)  = 

-.1288631602 

C  (  2 1)  = 

224.1213882024 

C ( 21) = 

70.9874042555 

C ( 22) = 

-355.1369302716 

C ( 22) = 

39.8818213389 

C ( 23) = 

115.5598816352 

C ( 23) = 

5.8251333617 

C  (  2  4)  = 

-82.6517911083 

C  ( 24)  = 

-10.3385864953 

C ( 25) = 

136.3420036065 

C  ( 25)  = 

-73.6285916503 

C (26) = 

.0000208297 

C ( 26) = 

.0000259077 

C (27) = 

3.6172294423 

C  ( 27)  = 

.0240818549 

C<28)  = 

-87.3706458740 

C ( 28 ) = 

-2.6011605493 

C  ( 29)  = 

69.9345875052 

C ( 29) = 

7.555445524 7 

C  ( 30)  = 

-10.0936503153 

C  (  30 )  = 

-.2349178007 

C  ( 31)  = 

-63.0570038645 

C  (  3 1 )  = 

-20.2504191554 

C  ( 32)  = 

105.8756257875 

C ( 32) = 

-13.2365012989 

C  ( 33)  = 

-35.1328556619 

C (33)= 

4.7647670677 

C  (  3  4 )  = 

28.2071799715 

C ( 34)= 

2.9103610175 

C ( 35) = 

-39  ?  68  5  0  2477 1 2 

C  ( 35)  = 

21.5041290069 

C  ( 36)  = 

-.0000087052 

C  (  36)  = 

-.0000121126 

C  ( 37)  - 

-.9340846773 

C  (  37)  = 

-.0138504245 

C  (  38)  = 

25.7701411555 

C ( 38) = 

-5.1011694772 

C ( 39) = 

-24.3306349045 

C  ( 39)  = 

-1.9502307115 

C(40)= 

4.5780233701 

C ( 40) = 

.1689371706 

i  - 96 -9 

66 %  within 

±  10% 

R!  =  97.8 
m 

65%  within  ± 

10% 

92%  within 

±  20% 

93%  within  ± 

20% 
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FIG.  32.  CRATER  VOLUME  AS  A  FUNCTION  OF  CHARGE  DEPTH 


AND  PERCENT  SATURATION,  S,  FOR  SPECIFIC  GRAVITY  =  2.55 


Scaled  Crater  Dimension,  L 


Dry  Unit  Weight,  y^,  Pounds/Cubic  Foot 

FIG,  36.  CRATER  RADIUS,  DEPTH  AND  VOLUME  AS  A  FUNCTION 
OF  DRY  UNIT  WEIGHT  FOR  OPTIMUM  CHARGE  DEPTH,  SPECIFIC 
GRAVITY  =2.55  AND  PERCENT  SATURATION  =  40 


1/3,,  ,5/1 6 


Scaled  Depth  of  Burst,  Dob/W^6 

FIG.  37.  CRATER  VOLUME  AS  A  FUNCTION  OF  CHARGE  DEPTH 
AND  DRY  UNIT  WEIGHT,  Yd,  POUNDS/CUBIC  FOOT,  FOR 
SPECIFIC  GRAVITY  =  2.65  AND  PERCENT  SATURATION  =  40 


145 


Percent  Saturation,  S 


FIG.  38.  CRATER  VOLUME  AS  A  FUNCTION  OF  PERCENT  SATURATION 
AND  DRY  UNIT  WEIGHT,  Yd,  POUNDS/CUBIC  FOOT,  FOR  OPTIMUM 
CHARGE  DEPTH  AND  SPECIFIC  GRAVITY  =  2.55 


REFERENCES 


1.  Baker,  Warren  J.,  “Effects  of  Soil  and  Rock  Properties 
on  Explosion  Crater  Characteristics,"  Master  of  Science 
dissertation,  Notre  Dame,  April  1962. 

2.  Banks,  D.C.,  and  Saucier,  R.T.,  "Geology  of  Buckboard 
Mesa,  Nevada  Test  Site,"  Report  PNE-5QQ1P,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg, 
Mississippi ,  July  1964. 

3.  Benfer,  Richard  H.,  Maj.,  "Project  Pre-Schooner  II, 
Apparent-Crater  Studies,"  Report  PNE-508,  U.S.  Army 
Engineer  Nuclear  Cratering  Group,  Livermore,  California, 

May  1967. 

4.  Benfer,  Richard  H.,  "Variation  of  Project  ZULU  II  Crater 
Dimensions  versus  Medium  Properties,"  NCG  Technical  Memo. 
66-18,  U.S.  Army  Engineer  Nuclear  Cratering  Group, 
Livermore,  California,  October  1966. 

5.  Bening,  Robert  G.,  Maj.,  and  Kurtz,  Maurice  K.  Jr.,  Lt, 
Col.,  "The  Formation  of  a  Crater  as  Observed  in  a  Series 
of  Laboratory-Scale  Cratering  Experiments,"  Report  PNE- 
5011 ,  U.S.  Army  Engineer  Nuclear  Cratering  Group, 

Livermore,  California,  October  1967. 

6.  Butkovich,  Theordore  R.,  "Calculation  of  the  Shock  Wave 
From  an  Underground  Nuclear  Explosion  in  Granite,"  Third 
Plowshare  Symposium  on  Engineering  with  Nuclear  Explosives, 
Report  TIP  7695,  Lawrence  Radiation  Laboratory,  LivermoreT 
California,  April  1964. 

7.  Carlson,  R.H.,  "Project  Toboggan,  High-Explosive  Ditching 
From  Linear  Charges,"  Research  Report  SC-4483  (RR),  Sandia 
Corp.,  Albuquerque,  New  Mexico,  July  1961. 

8.  Carlson,  R.H.,  (U)  "Review  of  Test  Results  and  Prediction 
Methods  Associated  with  Nuclear  Craters  and  Their  Related 
Ejecta  Distributions,"  Report  No.  D2-1 25706-1,  The  Boeing 
Co.,  Seattle,  Washington,  January  1968,  (Secret-RD). 

9.  Carlson,  R.H.,  and  Jones,  G.D.,  "Project  MTCE,  High 
Explosive  Cratering  Studies  in  Hard  Rock,"  Report  D2- 
90704-1,  The  Boeing  Co.,  Seattle,  Washington ,  November 
1965. 


147 


10.  Chabai ,  A.J.,  (U)  "Dimensional  Analysis  and  Modeling  of 
Cratering  Explosions,’1  Proceedings  of  the  Symposium  on 
Nuclear  Craters  and  Ejecta,  Report  tR-1 001  (52855-40) -4, 
Aerospace  Corp. ,  Novemoer  1967,  (U)  pp.  C-l-C-10,  (Secret- 
RD). 

11.  Chabai,  A.J.,  "Scaling  Dimensions  of  Craters  Produced  by 
Buried  Explosions,"  Research  Report  SC-RR-65-70,  Sandia 
Corp.,  Albuquerque,  New  Mexico,  February  1965. 

12.  Cherry,  J.T.,  "Computer  Calculations  of  Explosion-Produced 
Craters,"  International  Journal  of  Rock  Mechanics  and 
Mineral  Sciences  ,~Vol .  4,  1967,  pp.  1-22. 

13.  Circeo,  L.O.  Jr.,  "Nuclear  Excavation:  Review  and 
Analysis,"  Engineering  Geology ,  3,  Elsevier  Publishing 
Co.,  Amsterdam,  Netherl ands ,  June  1969,  pp.  5-59. 

14.  Cook,  Melvin  A.,  The  Science  of  High  Explosives,  Reinhold 
Publishing  Corp.,  New  York,  1958. 

15.  Dietz,  John  P.,  and  Hayes,  Dennis  B.,  "Compilation  of 
Crater  Data,"  Research  Report  SC-RR-65-220 ,  Sandia  Corp., 
Albuquerque,  New  Mexico,  July  1965. 

16.  Draper,  N.R.,  and  Smith,  H.,  Applied  Regression  Analysis, 
John  Wiley  and  Sons,  Inc.,  New  York,  1966 . 

17.  Fisher,  Paul  R.,  Kley,  Ronald  J.,  Capt.,  and  Jack,  Harold 
A.,  "Project  Pre-Gondola  I,  Geologic  Investigations  and 
Engineering  Properties  of  Craters,"  Report  PNE-1103,  U.S. 
Army  Engineer  Nuclear  Cratering  Group,  Livermore, 
California,  May  1969. 

18.  Flanagan,  T.J.,  "Project  Air  Vent,  Crater  Studies,"  Report 
SC-RR-64-1704 ,  Sandia  Corp.,  Albuquerque,  New  Mexico, 

April  1966. 

19.  Frandsen,  A. D. ,  "Analysis  and  Re-evaluation  of  Bulking 
Factors,"  TM  70-1 ,  U.S.  Army  Engineer  Nuclear  Cratering 
Group,  Lawrence  Radiation  Laboratory,  Livermore, 
California,  March  1970. 

20.  Frandsen,  A.D.,  "Engineering  Properties  Investigations  of 
the  Cabriolet  Crater,"  Report  PNE-957,  U.S.  Army  Engineer 
Nuclear  Cratering  Group,  Livermore,  California,  March  1970. 


148 


21.  Goode,  T.B.,  end  Mathews,  A.L.,  "Operation  Sun  Beam,  Shot 
Johnie  Boy,  Soils  Survey,  "  Report  PQR-2285,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg, 
Mississippi,  May  1963. 

22.  Hansen,  Spenst  M.,  et.  al.,  "Recommended  Crater  Nomen¬ 
clature,"  UCRL-7750 ,  University  of  California,  Lawrence 
Radiation  Laboratory ,  Livermore,  California,  March  1964. 

23.  Harlan,  R.W.,  Capt.,  "Project  Pre-Gondola  I,  Crater 
Studies:  Crater  Measurements,"  Report  PNE-1107,  U.S.  Army 
Engineer  Nuclear  Cratering  Group,  Livermore,  California, 
December  1967. 

24.  Hemmerle,  William  J.,  Statistical  Computations  on  a  Digital 
Computer,  Blaisdell  Publishing  Co.,  Waltham,  Massachusetts, 
1967. 

25.  Herr,  Robert  W.,  "A  Preliminary  Look  at  Effects  of 
Atmospheric  Pressure  and  Gravity  on  the  Size  of  Explosive 
Craters,"  Langley  Working  Paper-790,  NASA  Langley  Research 
Center,  Hampton,  Virginia,  August  1969. 

26.  Holzer,  Fred,  "Calculation  of  Seismic  Source  Mechanisms," 
Report  UCRL-1 221 9 ,  Lawrence  Radiation  Laboratory, 

Livermore,  California,  January  1965. 

27.  Hughes,  Bernard  C.,  Lt.  Col.,  "Nuclear  Construction 
Engineering  Technology,"  NCG  Technical  Report  No.  2,  U.S. 
Army  Engineer  Nuclear  Cratering  Group,  Livermore, 
California,  September  1968. 

28.  Hughes,  B.C.,  and  Benfer,  R.H.,  "Project  ZULU,  A  One- 
Pound  High  Explosive  Cratering  Experiment  in  Scalped  and 
Remoulded  Desert  Alluvium,"  NCG  Technical  Memo.  65-9, 

U.S.  Army  Engineer  Nuclear  Cratering  Group,  Livermore, 
California,  October  1965. 

29.  Hunt,  R.W.,  Bailey,  D.M.,  and  Carter,  L.D.,  "Preshot 
Geological  Engineering  Investigations  for  Project 
Cabriolet,  Pahute  Mesa,  Nevada  Test  Site,"  Report  PNE- 
966_,  U.S.  Army  Engineering  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  October  1968. 

30.  Johnson,  Gerald  W.,  et.  a'l . ,  "The  Plowshare  Program-History 
and  Goals,"  Third  Plowshare  Symposium  on  Engineering  with 
Nuclear  Explosives,  Report  Tip  7695.  Lawrence  Radiation 
Laboratory,  Livermore,  California,  ApriJ  1964. 


149 


31.  Johnson,  Gerald  W.,  et.  al.,  "The  Underground  Nuclear 
Detonation  of  September  19,  1957,  Ranier,  Operation 
Plumb-bob,"  Report  UCRL  5124,  University  of  California 
Radiation  Laboratory,  Livermore,  California,  February 
1958. 

32.  Knox,  J.B.,  and  Terhune,  R.W.,  "Calculation  of  Explosion- 
Produced  Craters,"  Third  Plowshare  Symposium  on  Engineering 
with  Nuclear  Explosives,  Report  TIP  7695,  Lawrence 
Radiation  Laboratory,  Livermore,  California,  April  1964. 

33.  Lutton,  R.J.,  et.  al . ,  "Project  Pre-Schooner  II,  Preshot 
Geologic  and  Engineering  Properties  Investigations," 

Report  PNE-509,  U.S.  Army  Engineer  Waterways  Exoeriment 
Station,  Vicksburg,  Mississippi,  December  1967. 

34.  Lutton,  R.J.,  and  Girucky,  F.E.,  "Project  Sulky,  Geologic 
and  Engineering  Properties  Investigations,"  Report  PNE- 
720,  U.S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  November  1966. 

35.  Lutton,  R.J.,  Girucky,  F.E.,  and  Hunt,  R.W.,  "Project 
Pre-Schooner,  Geologic  and  Engineering  Properties  Inves¬ 
tigations,"  Report  PNE-505F,  U.S.  Army  Engineer  Nuclear 
Cratering  Group,  Livermore,  California,  April  1967. 

36.  Maenchen,  George,  and  Nuckolls,  John,  "Calculation  of 
Underground  Explosions,"  Proceedings  of  the  Geophysical 
Laboratory  -  Lawrence  Radiation  Laboratory  Cratering 
Symposium,  Report  UCRL-6438,  Part  II,  Lawrence  Radi ati on 
Laboratory,  Livermore,  California,  October  1961,  pp. 
J-l-J-6. 

37.  Maenchen,  G.,  and  Sack,  S.,  "The  Tensor  Code,"  Report 
UCRL-7316 ,  Lawrence  Radiation  Laboratory,  Livermore, 
California,  April  1963. 

38.  "Military  Engineering  with  Nuclear  Explosives,"  Report 
DASA  1669,  U.S.  Army  Engineer  Nuclear  Cratering  Group, 
Livermore,  California,  June  1966. 

39.  Miller,  Irwin,  and  Freund,  John  E.,  Probability  and 
Statistics  for  Engineers,  Prentice-Hall ,  Inc.,  Englewood 
CTTffs,  New  Jersey,  1965. 

40.  Murphey,  Bryon  F.,  "High  Explosive  Crater  Studies:  Desert 
Alluvium,"  Research  Report  SC-4611  (RR),  Sandia  Corp., 
Albuquerque,  New  Mexico,  May  1961. 


150 


Murphey,  Bryon  F.,  "High  Explosive  Crater  Studies:  Tuff," 
Research  Report  SC-4574  (RR),  Sandia  Corp.,  Albuquerque, 
New  Mexico,  April  196T. 

Nelson,  David  L.,  1 Lt . ,  and  Taylor,  T.S.  Ill,  Sp4, 
"Analysis  of  Vesic  Crater  Modeling  Experiments,"  NCG 
Technical  Memo  68-14,  U.S.  Arn\y  Engineer  Nuclear  Cratering 
Group,  Livermore,  California,  November  1968. 

Nordyke,  M.D.,  "Nevada  Test  Site  Nuclear  Craters,"  Pro¬ 
ceedings  of  the  Geophysical  Laboratory  -  Lawrence 
Radiation  Laboratory  Cratering  Symposium,  Report  UCRL- 
6438,  Part  I,  Lawrence  Radiation  Laboratory,  Livermore, 
Cali  form  a,  "October  1961,  pp.  F-l-F-14. 

Nordyke,  M.D.,  Editor,  "Proceedings  of  the  Geophysical 
Laboratory  -  Lawrence  Radiation  Laboratory  Cratering, 
Symposium,"  Report  UCRL-6438,  Parts  I  &  II,  Lawrence 
Radiation  Laboratory ,  Livermore,  CaTTTornia,  October  1961. 

Nordyke,  M.D.,  "Cratering  Experience  with  Chemical  and 
Nuclear  Explosives,"  Third  Plowshare  Symposium  on 
Engineering  with  Nuclear  Explosives,  Report  TIP  7695, 
Lawrence  Radiation  Laboratory,  Livermore,  California, 

April  1964. 

Nordyke,  M.D.,  "On  Cratering,  A  Brief  History,  Analysis, 
and  Theory  of  Cratering,"  Report  UCRL-6578,  University  of 
California,  Lawrence  Radiation  Laboratory,  Livermore, 
California,  August  1961. 

Nordyke,  M.D.,  and  Williamson,  M.M.,  "Project  Sedan,  The 
Sedan  Event,"  Report  PNE-242F,  Lawrence  Radiation  Labo¬ 
ratory,  Livermore,  California,  August  1965. 

Nordyke,  M.D.,  and  Wray,  W.,  "Cratering  and  Radioactivity 
Results  from  a  Nuclear  Cratering  Detonation  in  Basalt," 
Journal  of  Geophysical  Research,  Vol .  69,  No.  4,  February 
15,  1964,  pp.  67'5-1689." 

Nugent,  R.C.,  and  Girucky,  F.E.,  "Project  Palanquin, 
Preshot  Geologic  and  Engineering  Properties  Investi¬ 
gations,"  Report  PNE-905,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  Mississippi,  October  1968. 


151 


50.  Nugent,  R.C.,  and  Banks,  D.C.,  "Project  Danny  Boy, 
Engineering  -  Geologic  Investigations,"  Report  PNE-5QQ5, 
U.S.  Army  Engineer  Nuclear  Cratering  Group,  Livermore, 
California,  November  1966. 

51.  Nugent,  R.C.,  and  Banks,  D.C.,  "Project  Pre-Schooner, 
Preshot  Investigations,"  Report  PNE-505P,  U.S.  Army 
Engineer  Nuclear  Cratering  Group,  Livermore,  California, 
September  1965. 

52.  Nugent,  R.C.,  and  Banks,  D.C.,  "Project  Sulky,  Preshot 
Geologic  Investigations,"  Report  PNE-719P,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg, 
Mississippi,  July  1965. 

53.  Perret,  William  R. ,  et.  a!.,  "Project  Scooter,"  Research 
Report  SC-4602  (RR) ,  Sandia  Corp.,  Albuquerque,  Riw 
Mexico,  October  1963. 

54 .  _ ,  Proceedings  of  the  Second  Plowshare  Symposium, 

Report  UCRL-5676,  Lawrence  Radiation  Laboratory,  Livermore, 
California,  and  AEC-San  Francisco  Operations  Office,  May 
1959. 

55.  _ ,  Proceedings  of  the  Symposium  on  Engineering  with 

Nuclear  Explosives,  Report  C0NF-700101,  Vols.  1  &  2,  The 
American  Nuclear  Society,  Las  Vegas,  Nevada,  May  1970. 

56.  Ramspott,  L.D.,  Stearns,  R.T.,  and  Meyer,  G.L.,  "Preshot 
Geophysical  Studies  of  Cabriolet  Nuclear  Crater  Site," 
Report  PNE-953  (UCRL-50458) ,  Lawrence  Radiation  Labo¬ 
ratory,  Livermore,  California,  August  1969. 

57.  Rooke,  A.D.  Jr.,  "Project  Pre-Buggy,  Emplacement  and 
Firing  of  High-Explosive  Charges  and  Crater  Measurements," 
Report  PNE-302,  U.S.  Army  Engineer  Nuclear  Cratering  Group, 
Livermore,  California,  February  1965. 

58.  Rooke,  A.D.  Jr.,  and  Davis,  L.R.,  "Ferris  Wheel  Series, 

Flat  Top  Event,  Crater  Measurements,"  Report  P0R-3008, 

U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 
Mississippi,  August  1966. 

59.  Sachs,  D.C.,  and  Swift,  L.M.,  "Small  Explosion  Tests, 
Project  Mole,"  Report  AFSWP-291 ,  Vols.  I  &  II,  Stanford 
Research  Institute,  Melo  Park,  California,  December  1955. 


152 


I 

t 

t 


60.  Sager,  R.A.,  Denzel,  C.W.,  and  Tiffany,  W.B.,  "Cratering 
from  High  Explosives,  Compendium  of  Crater  Data," 

Technical  Report  No.  2-547,  Report  1,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi, 

May  1960. 

61.  Saxe,  H.C.,  and  DelManzo,  D.D.  Jr.,  "A  Study  of  Underground 
Explosion  and  Cratering  Phenomena  in  Desert  Alluvium," 
Proceedings  of  the  Symposium  on  Engineering  with  Nuclear 
Explosives,  Report  C0NF-700l0r,  Vols.  T~1T~2,  The  American 
Nuclear  Society,  Las  Vegas,  Nevada,  May  1970,  pp.  1701- 
1725. 

62.  Saxe,  Harry  C.,  "Explosion  Crater  Prediction  Utilizing 
Characteristic  Parameters,"  Rock  Mechanics ,  C.  Fairhurst, 
ed.,  Pergaman  Press,  New  York,  19*6 "3. 

63.  Sauer,  F.M.,  Clark,  G.B.,  and  Anderson,  D.C.,  "Nuclear 
Geoplosives,  Part  Four-Empirical  Analysis  of  Ground  Motion 
and  Cratering,"  Report  DASA-1285  (IV),  Stanford  Research 
Institute,  Melo  ParlT,  California"  May  1964. 

64.  "Soil  Data,  Tonopah  Test  Site,  Nevada,"  Woodward-Clyde- 
Sherard  and  Associates,  report  to  Sandia  Corp.,  Denver, 
Colorado,  November  2.6,  1963. 

65.  Spruill,  Joseph  L.,  and  Paul,  Roger  A.,  "Project  Pre- 
Schooner,  Crater  Measurements,"  Report  PNE-502F,  U.S.  Army 
Engineer  Nuclear  Cratering  Group,  Livermore,  California, 

June  1965. 

66.  Spruill,  Joseph  L.,  and  Videon,  Fred  F.,  "Project  Pre- 
Buggy  II,  Studies  of  the  Pre-Buggy  II  Apparent  Craters," 
Report  PNE-315F,  U.S.  Army  Engineer  Nuclear  Cratering 
Group,  Livermore,  California,  June  1965. 

67.  Stearns,  R.T.,  Meyer,  G.L.?  and  Hansen,  S.M.,  "Project 
Palanquin,  Preshot  Geophysical  Properties  of  Palanquin 
Crater  Site,"  Report  PNE-906,  Lawrence  Radiation  Laboratory, 
Livermore,  California,  June  1968. 

68.  Stearns,  R.T.,  and  Rambo,  J.T.,  "Project  Pre-Gondola  I, 
Preshot  Geophysical  Measurements,"  Report  PNE-1111 ,  U.S. 

Army  Engineer  Nuclear  Cratering  Group,  Livermore, 

Cali fornia,  July  1967. 


153 


f 


i 


i 

i 


l 


69.  Strohm,  W.E.,  Ferguson,  J.S.  Jr.,  and  Krinitzsky,  E.L., 
"Project  Sedan,  Stability  of  Crater  Slopes,"  Report  PNE- 
234F,  U.S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  October  1964. 

70.  Teller,  Edward,  et.  al..  The  Constructive  uses  of  Nuclear 

Explosives ,  McGraw-Hill  .Book  Co.,  New  York,  1S68.  ~ 

71.  Terhune,  R.W.,  Stubbs,  T.F.,  and  Cherry,  J.T.,  "Nuclear 
Cratering  on  a  Digital  Computer,"  Report  UCRL-72032, 
Lawrence  Radiation  Laboratory,  Livermore,  California, 
February  1970. 

72.  Tewes,  Howard  A.,  "Results  of  the  Schooner  Excavation 
Experiment,"  Report  UCRL-71 832 ,  Lawrence  Radiation 
Laboratory,  Li vermore ,  Cali f orn i a ,  January  1970. 

73.  "The  Unified  Soil  Classification  System,"  Technical  Memo 
No.  3-357,  U.S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  March  1953,  (Rev.  April  1960). 

74.  Third  Plowshare  Symposium  on  Engineering  with  Nuclear 
Explosives,  Report  TIP  7695,  Lawrence  Radiation  Laboratory , 
Livermore,  California,  April  1964. 

75.  Thompson,  L.J.,  et.  al . ,  "Earth  Penetration  and  Dynamic 
Soil  Mechanics,  The  Effect  of  Soil  Parameters  on  Earth 
Penetration  of  Projectiles,"  Project  596-1,  report  to 
Sandia  Laboratories,  Texas  A&M  Research  Foundation,  College 
Station,  Texas,  July  1969. 

76.  "Underground  Explosion  Test  Program,  Sandstone,"  Technical 
Report  No.  5,  Vol .  I,  Engineering  Research  Associates, 

St.  Paul,  Minnesota,  February  1953. 

77.  Vaile,  R.B.  Jr. ,  "Pacific  Craters  and  Scaling  Laws," 
Proceedings  of  the  Geophysical  Laboratory  -  Lawrence 
Radiation  Laboratory  Cratering  Symposium,  Report  UCRL- 
6338,  Part  I,  Lawrence  Radiation  Laboratory,  Livermore, 
California,  October  1961,  pp.  E-l-E-36. 

78.  Vesic,  A.S.,  "Cratering  by  Explosives  as  an  Earth  Pressure 
Problem,"  Sixth  International  Conference  on  Soil  Mechanic 
and  Foundation  Engineering,  Vo).  Il ,  "Montreal ,  Canada, 

1965,  pp.  427-431. 


154 


X 

i 


l 


<3 


J 


79.  Vesic,  A.S.,  Soutwell,  G.P.  Jr.,  and  Tai ,  Tein-Lie, 
"Engineering  Properties  of  Nuclear  Craters,  Theoretical 
Studies  of  Cratering  Mechanisms  Affecting  the  Stability 

of  Cratered  Slopes,  Phase  III,"  Technical  Report  No.  3-699, 
Report  6,  U.S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  March  1967. 

80.  Vesic,  A.S.,  et.  al . ,  "Engineering  Properties  of  Nuclear 
Craters,  Theoretical  Studies  of  Cratering  Mechanisms 
Affecting  the  Stability  of  Cratered  Slopes,  Phase  II," 
Technical  Report  No.  3-699,  Report  2,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi, 
October  1965. 

81.  Vesic,  Aleksandar  B.,  and  Barksdale,  Richard  D. , 
"Theoretical  Studies  of  Cratering  Mechanics  Affecting 
the  Stability  of  Cratered  Slopes,"  Final  report.  Project 
No.  A-655,  Engineering  Experiment  Station,  Georgia 
Institute  of  Technology,  Atlanta,  Georgia,  September  1963. 

82.  Videon,  F.F.,  "Project  Palanquin,  Studies  of  Apparent 
Crater,"  Report  PNE-904,  U.S.  Army  Engineer  Nuclear 
Cratering  Group,  Livermore,  California,  July  1966. 

83.  Videon,  F.F.,  "Project  Sulky,  Crater  Measurements,"  Report 
PNE-713F,  U.S.  Army  Engineer  Nuclear  Cratering  Group, 
Livermore,  California,  October  1965. 

84.  Vortman,  L.J.,  "Craters  from  Surface  Explosions  and 
Scaling  Laws"  (Reprint),  Sandia  Laboratories,  Albuquerque, 
New  Mexico,  July  1968. 

85.  Vortman,  L.J.,  et.  al.,  "Project  Buckboard,  20-Ton  and 
1/2-Ton  High  Explosive  Cratering  Experiments  in  Basalt 
Rock,"  Report  SC-4675  (RR),  Sandia  Corp.,  Albuquerque, 

New  Mexico,  August  1963. 

86.  Vortman,  L.J.,  et.  al.,  "Project  Stagecoach,  20-Ton  HE 
Cratering  Experiments  in  Desert  Alluvium,"  Report  SC- 
4596  (RR),  Sandia  Corp.,  Albuquerque,  New  Mexico,  March 

m: 

87.  Vortman,  L.J.,  "Ten  Years  of  High  Explosive  Cratering 
Research  at  Sandia  Laboratory,"  Nuclear  Applications  & 
Technology,  Vol .  7,  September  1969,  pp.  269-304. 

88.  Westine,  P.S.,  "Explosive  Cratering,"  Journal  of  Terra- 
mechanics  ,  Vol.  7,  No.  2,  Pergamon  Press,  1970,  pp.  9-19. 


155 


89.  Whitman,  Robert  V.,  "Soil  Mechanics  Consideration 
Pertinent  to  Predicting  the  Immediate  and  Eventual  Size 
of  Explosion  Craters  ,"  Report  SC-4404  (RR),  Sandia  Corp. , 
Albuquerque,  New  Mexico,  December  1959. 

90.  Whitman,  Robert  V.,  and  Healy,  Kent  A.,  "The  Response  of 
Soils  to  Dynamic  Loadings,  Report  9:  Shearing  Resistance 
of  Sands  During  Rapid  Loadings,"  Research  Report  R-62-113 
Massachusetts  Institute  of  Technology,  Cambridge,  “ 
Massachusetts,  May  1962. 

\ 

91.  Wu,  Tein  Hsing,  Soil  Mechanics,  Allyn  and  Bacon,  Inc., 

Boston,  1966.  / 


UNCLASSIFIED _ 

^  .  ?>curitj'_Clas«ific»1iori 

DOCUMENT  CONTROL  DATA  ■  R  4  D 

. (Security,  classification  ot  title  body  of  abstract  and  indexing  annotation  must  be  entered  uhen  the  overall  report  Js  classified) 

I*  ORIGIN  A  TINC  ACTIVITY  (Corporate  author)  (2*.  REPOrV  SC(i  JRITV  CLASSIFICATION 

•  Texas  AJ.M  Research  Foundation  UNCLASSIFIED 

College  Station,  Texas  77843  |26,  CROUP 


oM-T^II  TL  Ej  ..  . 


JHE  INFLUENCE  OF  SOIL  AND  BOCK  PROPERTIES  ON  THE  DIMENSIONS  OF  / 
'EXPLOSION-PRODUCED  CRATERS?  ~  / 

.  —  7 ffr-3 — ; — rzrizzzL - - — . — 

I  4.  DESCRIPTIVE  NOTES  l,7Vo»  ol  report  Mild  inclusive  dmleM)  /  t  _  / 

Fe<mKH>-y-4$7fln^  H  UcLw(„.  I  A  Jlit  h+b  '?$- Qci 


f  ^rjj^THORlSi  (First  name,  middlo  Initial,  last  niint*; 

tarry  A./DiTTonl  Maj,  USAF 


•  a .  CONTRACT  OR  GRAN] 


If/ 


CmMki  _ 

la.  TOTAL  PACE5  NO.  OF 

DA^QRIGINyOIVS  R  EPQ  RT^UM  PL  3(S) 

(  V^AFWLWE-71-144  i 

,  J^ifczrz _ 

[  05.  OTHER  REPORT  NOtSl  (Any  other  numbers  that  may  b»  resigned 

I  this  report) 


10.  DISTRIBUTION  STATEMENT 


Distribution  limited  to  US  Government  agencies  only  because  of  test  and  evaluation 
(19  Jan  72).  Other  requests  for  this  document  must  be  referred  to  AFWL  (DEV). 


II.  SUPPLEMENTARY  NOTES 


12*  SPONSORING  MILITARY  ACTIVITY 

AFWL  (DEV) 

Kirtland  AFB,  NM  87117 


[is.  abstract  ( D i s t r 1  b u t i on  Limitation  Statement  B) 

Analysis  of  data  from  published  cratering  experiments  shows  the  effect  of  soil  and 
rock  properties  on  the  apparent  dimensions  of  explosion-produced  craters.  More  than 
200  cratering  tests  and  related  material  properties  were  cataloged.  The  data  con¬ 
sisted  of  10  nuclear  events  whose  yields  varied  from  0.42  to  100  kilotons  and  about 
200  high  explosive  events  whose  yields  varied  from  1  to  1  million  pounds  of  TNT.  The 
different  test  sites  included  materials  for  which  the  density  ranged  from  60  to  170 
pounds/cubic  foot.  By  regression  analysis,  using  bell  shaped  curves,  prediction 
formulas  were  developed  for  the  apparent  crater  radius,  depth,  and  volume  as  a  func¬ 
tion  of  charge  weight  and  depth  of  burst  for  eight  different  types  of  materials.  The 
bell  curves  were  normalized  using  material  properties  and  prediction  equations  were 
generated  using  all  the  data.  These  general  equations  were  then  studied  to  determine 
the  specific  effects  of  the  material  properties  on  resultant  apparent  crater  dimen¬ 
sions.  Material  properties  are  highly  important  in  determining  the  size  of  explo¬ 
sion-produced  craters  and  some  of  the  more  important  properties  are  unit  weight, 
degree  of  saturation,  shearing  resistance  and  seismic  velocity.  Previous  investiga¬ 
tors  have  been  somewhat  negligent  in  measuring  material  properties  for  past  cratering 
experiments  and  no  real  data  analysis  can  be  made  until  the  variables  are  either 
controlled  or  measured.  Material  properties  which  should  be  measured  for  future 
tests  should  at  least  include  the  above  properties  and  if  possible  the  material's 
energy  dissipation  and  bulking  characteristics.  Better  yet  a  reasonably  simple 
theory  of  cratering  is  needed  which  will  better  define  the  material  properties 
governing  cratering  mechanics. 
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